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Abstract: We study the dispersion relations of mesons in a particular hot strongly coupled su- 
persymmetric gauge theory plasma. We find that at large momentum k the dispersion relations 



become lu ~ t> /c + a + b/k + . . ., where the limiting velocity vq is the same for mesons with any 
quantum numbers and depends only on the ratio of the temperature to the quark mass T/m q . We 
compute a and b in terms of the meson quantum numbers and T/m q . The limiting meson veloc- 
t^. ' ity vq becomes much smaller than the speed of light at temperatures below but close to Tdi SS , the 
temperature above which no meson bound states at rest in the plasma are found. From our result 
for Vo(T/m q ), we find that the temperature above which no meson bound states with velocity v 
exist is T diss (t>) ~ (1 — t ,2 ) 1//4 T diss , up to few percent corrections. We thus confirm by direct calcu- 
lation of meson dispersion relations a result inferred indirectly in previous work via analysis of the 
screening length between a static quark and antiquark in a moving plasma. Although we do not 
do our calculations in QCD, we argue that the qualitative features of the dispersion relation we 
compute, including in particular the relation between dissociation temperature and meson velocity, 
may apply to bottomonium and charmonium mesons propagating in the strongly coupled plasma 
of QCD. We discuss how our results can contribute to understanding quarkonium physics in heavy 
ion collisions. 
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1. Introduction 



The radii of the tightly bound heavy quark- ant iquark systems of the charmonium (J/^f, Xc, •••) 
and bottomonium (T, T', ...) families provide a unique set of decreasing length scales in strong 
interaction physics. On general grounds, it is expected that the attraction between a heavy quark 
and an anti-quark is sensitive to the medium in which the bound state is embedded, and that 
this attraction weakens with increasing temperature. In the context of ultra-relativistic nucleus- 
nucleus collisions, the radii of some quarkonia states correspond to fractions of the natural length 
scale displayed by the medium produced in heavy ion collisions, namely fractions of its inverse 
temperature 1/T. Such scale considerations support the idea that measurements of the medium- 
modification or dissociation of quarkonia can characterize properties of the QCD matter produced 
in heavy ion collisions. 

Matsui and Satz were the first to highlight the role of quarkonium in the study of hot QCD 
matter |TJ. They suggested that J/\l/-suppression is a signature for the formation of deconfined 
quark-gluon plasma (QGP). More precisely, they argued that in comparison to proton-proton or 
proton-nucleus collisions, the production of J/^f mesons should be suppressed if quark-gluon plasma 
is formed in sufficiently energetic nucleus-nucleus collisions, since the screened interaction of a c 
and a c in QGP would not bind them [|IJ. The theoretical basis for this argument has been clarified 
considerably within the last two decades . Model-independent calculations of the static potential 
between a heavy quark and anti-quark have been performed in lattice-regularized QCD, valid at 
strong coupling |3|, |j, ^ || [7|, In lattice calculations without dynamical quarks, at temperature 
T = and large separation L this potential rises linearly with L, consistent with confinement. 
At nonzero temperature, the potential weakens and levels off at large distances; with increasing 
temperature, the distance at which this screening occurs decreases. This behavior of the static po- 
tential has been mapped out for hot QCD matter both without [|J and with || ||] dynamical quarks. 
However, the physical interpretation of static potentials at finite temperature rests on additional 
assumptions. For instance, even if a potential supports a bound state with several MeV binding 
energy, it remains unclear which physics can be attributed to such a state in a heat bath of ~ 200 
MeV temperature. Such issues do not arise in a discussion of quarkonium mesons based directly on 
their Minkowski space spectral functions or dispersion relations. In recent years, the spectral func- 
tions have been characterized by lattice calculations of the Euclidean correlation functions to which 
they are analytically related, again in hot QCD matter both without || and with [IIJ dynamical 
quarks. The use of these calculations of finitely many points on a Euclidean correlator to con- 
strain the Minkowski space spectral function of interest via the Maximum Entropy Method requires 
further inputs — for example smoothness assumptions or information on the analytic properties 
of the spectral function || jlQfl . At high enough temperatures that quark-gluon plasma becomes 
weakly coupled, a complementary analytical approach based upon resummed hard-thermal-loop 
perturbation theory becomes available ||11|| . These calculations have the advantage that analyti- 
cal continuation from Euclidean to Minkowski space does not introduce additional uncertainties, 
but it remains unclear to what extent they can treat a strongly coupled quark-gluon plasma. In 
broad terms, all these calculations support the qualitative picture behind the original suggestion 
of Matsui and Satz that color screening in the quark-gluon plasma is an efficient mechanism for 
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quarkonium dissociation. In addition, these studies support the picture of a sequential dissociation 
pattern in which loosely bound, large, quarkonia such as the and Xc cease to exist close 



to T c , the temperature of the crossover between hadronic matter and quark-gluon plasma, whereas 
more tightly bound, smaller, states dissociate only at significantly higher temperatures. In partic- 
ular, J/ty mesons continue to exist for a range of temperatures above the QCD phase transition 



and dissociate only above a temperature that lies between 1.5 T c and 2.5 T c |12| . The observation of 
bound-state-specific quarkonia suppression patterns could thus provide detailed information about 
the temperature attained in heavy ion collisions. 

On the experimental side, there are by now data from the NA50 and NA60 experiments at the 
CERN SPS and from the PHENIX experiment at RHIC demonstrating that the production of J/^> 
mesons is suppressed in ultra-relativistic nucleus-nucleus collisions compared to proton-proton or 
proton-nucleus collisions at the same center of mass energy [O. However, due to lack of statistics 



and resolution, an experimental characterization of other charmonium states (^f', Xc, • ••) has not 
yet been possible at RHIC, and bottomonium states have not yet been characterized in any nucleus- 
nucleus collisions. Moreover, the observed yield of J/\l/ mesons is expected to receive significant 
decay contributions from ty' and Xc meaning that the observed suppression of J/\l/ mesons may 



originate only in the suppression of the larger and Xc states f[2" |, or may indicate a suppression 
in the number of primary J/\l/ mesons themselves in addition. Thus, at present an experimental 
test of the sequential quarkonium suppression pattern is not in hand. It is expected that the LHC 



heavy ion program will furnish such a test, since two LHC experiments |14| have demonstrated 



capabilities for discriminating between the different states of the charmonium and bottomonium 
families. 

From the existing data in ultra-relativistic heavy ion collisions and their phenomenological 
interpretation, it has become clear that an unambiguous characterization of color screening effects in 
the quarkonium systems requires good experimental and theoretical control of several confounding 
factors. These include in particular control over the spatio-temporal evolution of the medium, 
control over the time scale and mechanism of quarkonium formation, as well as control over the 
effects of quarkonium propagation through the medium. We now comment on these three sources 
of uncertainty in more detail: 

First, there is ample evidence by now that the systems produced in ultra-relativistic heavy 
ion collisions display effects of position-momentum correlated motion (a.k.a. flow), which are as 
important as the effects of random thermal motion [[TJJ. Moreover, the energy density achieved in 
these collisions drops rapidly with time as the matter expands and falls apart after approximately 
10 fm/c. As a consequence, the modeling of quarkonium formation in heavy ion collisions cannot 
be limited to a description of heavy quark bound states in a heat bath at constant temperature 
(which is the information accessible in ab initio calculations in lattice-regularized QCD). The effects 
of a rapid dynamical evolution during which the relevant degrees of freedom in the medium change 
from partonic to hadronic must be taken into account. 

Second, regarding the formation process, the conversion of a heavy quark pair produced in 
a hard collision into a bound quarkonium state is not fully understood, even in the absence of a 
medium. There are different production models, which all have known limitations and for which 
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a systematic calculation scheme remains to be fully established (for a short review of these issues, 
TBI). The need for further clarification of the vacuum case has even led to suggestions that nu- 



scc 



clear matter could serve as a filter to distinguish between different production mechanisms |fl7 
However, it has also been pointed out that there may be a novel quarkonium production mechanism 
operating only in ultra- relativistic heavy ion collisions at RHIC and at the LHC charm quarks 
may be so abundant in these collisions that c and c quarks produced separately in different primary 
hard scattering interactions may find each other and combine, contributing significantly to charmo- 
nium production at soft and intermediate transverse momentum. To a lesser extent, this mechanism 
may also contribute to the production of Upsilon mesons. Identifying and characterizing such a 
novel formation process is of considerable interest, since recombination is likely to be quadratically 
sensitive to the phase space density of charm and thus to properties of the produced matter. On 
the other hand, if realized in nature recombination also implies that quarkonium spectra at soft 
and intermediate transverse momenta are determined predominantly during the late hadronization 
stage and cannot be viewed as probes which test color screening in the quark gluon plasma. This 
would indicate that the high transverse momentum regime (say above 5-8 GeV) of quarkonium 
spectra, which should not be significantly affected by recombination, is better suited for tests of the 
fundamental color screening effects predicted by QCD. However, the sensitivity of high transverse 
momentum spectra to properties of the medium remains to be established. In particular, quarko- 
nium formation or dissociation proceeds on a time scale comparable to the size of the bound state 
in its rest frame, meaning that quark-antiquark pairs with very high transverse velocity may escape 
the finite-sized droplet of hot matter produced in a heavy ion collision before they have time to 
form a meson, meaning in turn that screening effects cease to play a role in quarkonium production 



above some very high transverse momentum [20|. At lower transverse momenta, where screening 



does play a role, one must nevertheless understand for how long quarkonium is exposed to the 
medium and how readily it dissociates if moving relative to that medium. For quarkonium at high 
transverse momentum, the time of exposure to the medium depends on the geometry of the collision 
region, which determines the in-medium path length, and it depends on the propagation velocity. 
The results contained in this paper give novel input to modeling this process by demonstrating that 
the real part of the finite temperature quarkonium dispersion relation can differ significantly from 
the vacuum one, and can imply a limiting quarkonium propagation velocity which is much smaller 
than c, the velocity of light in vacuum. Our results indicate that at temperatures close to but below 
that at which a given quarkonium state dissociates, these mesons move through a strongly coupled 
quark-gluon plasma at a velocity that is much smaller than c even if they have arbitrarily high 
transverse momentum. Certainly this means that the formation time arguments of |2fJ will need 
rethinking before they can be applied quantitatively. 

Third, we turn to the question of how the relative motion of quarkonium with respect to the 
local rest frame of the medium affects quarkonium production. As discussed above, the standard 
vacuum relation between the momentum of a quarkonium state and its velocity can be altered in the 
presence of a medium and this effect may be phenomeno logically relevant. In addition, it is expected 
that a finite relative velocity between the medium and the bound state enhances the probability of 
dissociation [21]. In a recent strong coupling calculation of hot M = 4 supersymmetric QCD, three 



of us have have shown [22, 23] that the sceening length L s for a heavy quark-antiquark pair decreases 
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with increasing velocity as L s (v,T) ~ L s (0, T)/^Fj, with 7 = l/yl — ^ 2 the Lorentz boost factor. 
This suggests that a quarkonium state that is bound at v = at a given temperature could dissociate 
above some transverse momentum due to the increased screening, providing a significant additional 
source of quarkonium suppression at finite transverse momentum. The present work started from 
the motivation to establish how this velocity scaling manifests itself in a description of mesons at 
finite temperature, rather than via drawing inferences from a calculation of the screening length 
that characterizes the quark-antiquark potential. This motivation is analogous to that behind going 
from lattice QCD calculations of the static potential in QCD to calculations of the Minkowski space 
meson spectral function. We shall do our calculation in a different strongly coupled gauge theory 
plasma, in which we are able to do this investigation for mesons with nonzero velocity. We shall see 
that the critical velocity for the dissociation of quarkonium inferred from the velocity scaling of the 
screening length also appears as a limiting velocity for high-momentum quarkonium propagation in 
the hot non-abelian plasma. 

Finally, the characterization of color screening also depends on the experimental and theoretical 
ability to separate its effects on quarkonium production from effects arising during the late time 
hadronic phase of the heavy ion collision. In particular, it has been noted early on that significant 
charmonium suppression may also occur in confined hadronic matter |23j. However, it has been 
argued on the basis of model estimates for the hadronic J/ty dissociation cross section [^] that 
dissociation in a hadronic heat bath is much less efficient than in a partonic one. The operational 
procedure for separating such hadronic phase effects is to measure them separately in proton-nucleus 



collisions ||26|| , and to establish then to what extent the number of J /\& mesons produced in nucleus- 
nucleus collisions drops below the yield extrapolated from proton- nucleus collisions fl3| , |27j . 

The above discussion highlights the extent to which an understanding of quarkonium pro- 
duction in heavy ion collisions relies on theoretical modelling as the bridge between experimental 
observations and the underlying properties of hot QCD matter. This task involves multiple steps. 
It is of obvious interest to validate or constrain by first principle calculations as many steps as 
possible, even in a simplified theoretical setting. The present work is one of a number of recent 
developments [[28|] that explore to what extent techniques from string theory, in particular the 
AdS/CFT correspondence, can contribute to understanding processes in hot QCD by specifying 
how these processes manifest themselves in a large class of hot strongly coupled non-abelian gauge 
theories. Although it is not known how to extend the AdS/CFT correspondence to QCD, there 
are several motivations for turning to this technique. First, there are a growing number of explicit 
examples which indicate that a large class of thermal non-abelian field theories with gravity duals 
share commonalities such that their properties in the thermal sector are either universal at strong 
coupling, i.e. independent of the microscopic dynamics encoded in the particular quantum field 
theory under study, or their properties are related to each other by simple scaling laws e.g. de- 
pending on the number of elementary degrees of freedom. This supports the working hypothesis 
that by learning something about a large class of strongly coupled thermal non-abelian quantum 
field theories, one can gain guidance towards understanding the thermal sector of QCD. Second, 
the AdS / CFT correspondence allows for a technically rather simple formulation of problems involv- 
ing real-time dynamics. This is very difficult in finite temperature lattice-regularized calculations, 
which exploit the imaginary time formalism. In particular, this is the reason why so far lattice 
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QCD calculations treat only static quark-antiquark pairs in the plasma, and why the only nonper- 
turbative calculation of the velocity dependence of quarkonium dissociation exploits the AdS/CFT 
correspondence. Third, data from experiments at RHIC pertaining to many aspects of the matter 
produced in heavy ion collisions indicate that this matter is strongly coupled. Since the AdS / CFT 
correspondence provides a mapping of difficult nonperturbative calculations in a quantum field 
theory with strong coupling onto relatively simple, semi-classical calculations in a gravity dual, it 
constitutes a novel — and often the only — technique for addressing dynamical questions about 
hot strongly coupled non-abelian matter, questions that are being raised directly by experimental 
results on QCD matter coming from RHIC. 

We have focussed in this Section on the larger context for our results. In Section 2, which is 
an introduction in a more narrow sense, we review the past results which serve as an immediate 
motivation for our work, in particular the screening length that characterizes the potential between 
a static quark and antiquark in a moving plasma wind. Adding fundamental quarks with finite mass 
m q , and hence mesons, into M = 4 SYM theory requires adding a D7-brane in the dual gravity 
theory, as we review in Section 3. The fluctuations of the D7-brane are the mesons, as we review 
for the case of zero temperature in Section 3. In Section 4 we set up the analysis of the mesons at 
nonzero temperature, casting the action for the D7-brane fluctuations in a particularly geometric 
form, written entirely in terms of curvature invariants. Parts of the derivation are explained in 
more detail in Appendix A. With all the groundwork in place, in Section 5 we derive the meson 
dispersion relations. In addition to obtaining them numerically without taking any limits as has 
been done previously pSfl , we are able to calculate them analytically in three limits: first, upon 
taking the low temperature limit at fixed k; second, upon taking the low temperature limit at fixed 
kT; and third, using insights from the first two calculations, at large k for any temperature. At 
large k we find 



where vq is independent of meson quantum numbers, depending only on T/m q . Vq turns out to 
be given by the local speed of light at the "tip of the D7-brane", namely the place in the higher 
dimensional gravity dual theory where the D7-brane comes closest to the black hole ]29|]. We 
compute a and b in terms of meson quantum numbers and T/m q . Our result for the limiting velocity 
Vq for mesons at a given temperature T can be inverted, obtaining T diss (u), the temperature above 
which no mesons with velocity v are found. We find that up to few percent corrections, our result 
can be summarized by 



where Tdi SS is the temperature at which zero-velocity mesons dissociate, obtained in previous work 
and introduced in Section 3. As we discuss in Section 2, our results obtained by direct calculation 
of meson dispersion relations confirm inferences reached (in two different ways) from the analysis 
of the screened potential between a static quark and antiquark in a hot plasma wind. In Section 
6, we close with a discussion of potential implications of these dispersion relations for quarkonia 
in QCD as well as a look at open questions. The dispersion relations that we calculate in this 
paper describe how mesons propagate and so affect a class of observables, but determining whether 
quarkonium meson formation from a precursor quark-antiquark pair is suppressed by screening is a 



lo = v k + a + — + . . . 
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more dynamical question that can at present be addressed only by combining our calculation and 
the more heuristic results of [E3] . 



2. From screening in a hot wind to moving mesons 

In the present work, we shall use the AdS / CFT correspondence to study the propagation of mesonic 
excitations moving through a strongly coupled hot quark-gluon plasma. In this Section, however, 
we introduce what we have learned from the simpler calculation of the potential between a test 
quark-antiquark pair moving through such a medium. This will allow us to pose the questions that 
we shall address in the present paper. 

The simplest example of the AdS/CFT correspondence is provided by the duality between 
Af = 4 super Yang-Mills (SYM) theory and classical gravity in AdS$ x 5*5 f30fl . A/" = 4 super 
Yang-Mills (SYM) theory is a conformally invariant theory with two parameters: the rank of the 
gauge group N c and the 't Hooft coupling A = g\ M N c . In the large N c and large A limit, gauge 
theory problems can be solved using classical gravity in AdS$ x S5 geometry. We shall work in this 
limit throughout this paper. 

In M = 4 SYM theory at zero temperature, the static potential between a heavy external quark 
and antiquark separated by a distance L is given in the large N c and large A limit by [BT, 



where the 1/L behavior is required by conformal invariance. This potential is obtained by computing 
the action of an extremal string world sheet, bounded at r — > 00 (r being the fifth dimension of 
AdSs) by the world lines of the quark and antiquark and "hanging down" from these world lines 



toward smaller r. At nonzero temperature, the potential becomes |33 

V(L,T) « VXf(L) L<L C 

« \°g(L) L> L c . (2.2) 

In ( |2.2| ), at L c = 0.24/T there is a change of dominance between different saddle points and the 
slope of the potential changes discontinuously. When L < L c , the potential is determined as 
at zero temperature by the area of a string world sheet bounded by the worldlines of the quark 
and antiquark, but now the world sheet hangs down into a different five-dimensional spacetime: 
introducing nonzero temperature in the gauge theory is dual to introducing a black hole horizon 
in the five-dimensional spacetime. When L <C L c , f(L) reduces to its zero temperature behavior 
( En]) . When L > L c , g(L) has the behavior [|4 



g{L) oc ci - c 2 e- m ^ L , (2.3) 

with ci, C2 and m gap constants all of which are proportional to T. This large-L potential arises 
from two disjoint strings, each separately extending downward from the quark or antiquark all the 
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way to the black hole horizon, exchanging supergravity modes the lightest of which has a mass 
given by m gap = 2.34 ttT. (There are somewhat lighter modes with nonzero R-charge, but these are 



not relevant here [35].) It is physically intuitive to interpret L c as the screening length L s of the 
plasma since at L c the qualitative behavior of the potential changes. Similar criteria are used in 
the definition of screening length in QCD ||, although in QCD there is no sharply defined length 
scale at which screening sets in. Lattice calculations of the static potential between a heavy quark 
and antiquark in QCD indicate a screening length L s ~ 0.5/T in hot QCD with two flavors of light 
quarks || and L s ~ 0.7/T in hot QCD with no dynamical quarks [|J. The fact that there is a 
sharply defined L c in fl2.2p is an artifact of the limit in which we are working. 1 

In |22] , |23|| , three of us studied the velocity scaling of the screening length L s in M = 4 super- 
Yang-Mills theory and found that 2 

L s (vAT) = ^(l-v>) 1/ \ (2.4) 

where 9 is the angle between the orientation of the quark- antiquark dipole and the velocity of the 
moving thermal medium in the rest frame of the dipole. f(v, 6) is only weakly dependent on both 
of its arguments. That is, it is close to constant. So, to a good approximation we can write 

L s (v, T) « L.(0, T) (1 - v 2 ) V* oc 1(1 - v 2 ) 1 " . (2.5) 



This result, also obtained in [|36[ and further explored in |37|, |39| , has proved robust in the 



sense that it applies in various strongly coupled plasmas other than Af = 4 SYM |37|, |38], |39| . The 
velocity dependence of the screening length (|2.5|) suggests that in a theory containing dynamical 
heavy quarks and meson bound states (which M = 4 SYM does not) the dissociation temperature 
Tdiss(v), defined as the temperature above which mesons with a given velocity do not exist, should 
scale with velocity as 

T diss (v) ~T diss (v = 0)(l-v 2 ) 1/i , (2.6) 
since T diss (v) should be the temperature at which the screening length L s (v) is comparable to the 



size of the meson bound state. The scaling (2J3) then indicates that slower mesons can exist up to 
higher temperatures than faster ones. In this paper, we shall replace the inference that takes us 
from the calculated result (|2.5|) to the conclusion (|2.6|) by a calculation of the properties of mesons 



1 The theoretical advantage of using l/m gap to define a screening length as advocated in |34| is that it can be 
precisely defined in JV = 4 SYM theory at finite A and N c , as well as in QCD, as it characterizes the behavior of the 
static potential in the L —> oo limit. The disadvantage of this proposal from a phenomenological point of view is that 
quarkonia are not sensitive to the potential at distances much larger than their size. For questions relevant to the 
stability of bound states, therefore, the length scale determined by the static potential that is phenomenologically 
most important is that at which the potential flattens. Although this length is not defined sharply in QCD, it is 
apparent in lattice calculations and can be defined operationally for practical purposes [W o. This L s seems most 



analogous to L c in ( 2J2 ) , and we shall therefore continue to refer to L a = L c as the screening length, as in the original 
literature [p3| . Note that L c is larger than l/m gap by a purely numerical factor ~ 1.8. 

2 In |2^, |23| L s was defined using a slightly different quantity than L c in (^2), such that L s = 0.28/T for a 



quark-antiquark at rest. For technical reasons, this other definition was more easily gcncralizable to nonzero velocity. 
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themselves, specifically their dispersion relations. We shall reproduce (|2.6|) in this more nuanced 
setting, finding few percent corrections to the basic scaling result inferred previously. 

The results fl2.5| ) and ( |2.6| ) have a simple physical interpretation which suggests that they 



could be applicable to a wide class of theories regardless of specific details. First, note that since 
L s (0) ~ ~, both ( |2.5| ) and Q2.6| ) can be interpreted as if in their rest frame the quark- ant iquark 
dipole experiences a higher effective temperature Although this is not literally the case 

in a weakly coupled theory, in which the dipole will see a redshifted momentum distribution of 
quasiparticles coming at it from some directions and a blueshifted distribution from others [^l|], we 
give an argument below for how this interpretation can nevertheless be sensible. The result ( |2.5| ) 
can then be seen as validating the relevance of this interpretation in a strongly coupled plasma. The 
argument is based on the idea that quarkonium propagation and dissociation are mainly sensitive 
to the local energy density of the medium. Now, in the rest frame of the dipole, the energy density 
(which we shall denote by p) is blue shifted by a factor ~ 7 2 and since p oc T 4 in a conformal 
theory, the result ( |2.5|) is as if quarks feel a higher effective temperature given by T^/7. 3 Lattice 
calculations indicate that the quark-gluon plasma in QCD is nearly conformal over a range of 
temperatures 1.5T C < T < 5T C , with an energy density p « bT A where b is a constant about 80% of 



the free theory value P0[ . So it does not seem far-fetched to imagine that (|2.5|) could apply to QCD. 
We should also note that AdS/CFT calculations in other strongly coupled gauge theories with a 
gravity description are consistent with the interpretation above [38[] and that for near conformal 



theories the deviation from conformal theory behavior appears to be small |[38|| . If a velocity scaling 
like ( |2.5| ) and ( |2.6| ) holds for QCD, it can potentially have important implications for quarkonium 
suppression in heavy ion collisions, as we have discussed in Section 1 and will return to in Section 6. 

While the argument leading from ( [2.5| ) to ( |2.6| ) is plausible, it is more satisfying to have a 
set-up within which one can study mesons directly. Direct study of meson bound states will also 
yield more insights than the study of the screening length from the potential. It is the purpose of 
this paper to examine this issue in a specific model with dynamical flavors. 

Before beginning our analysis, let us first note a curious feature regarding the quark potential 
observed in [22, 23| . There one introduces a probe brane near the boundary of the AdSs black 



hole geometry with open strings ending on it corresponding to fundamental "test quarks" of mass 
given by 



V\T. It was found that for any given quark mass m„, there exists a maximal velocity v c 



= 1 - si} • (27) 

beyond which there is no 0{\f\) potential between the pair for any value of their separation larger 
than their Compton wavelength, i.e. for any distance at which a potential can be defined. This 
result can be interpreted as saying that for any given T and m q , it is impossible to obtain bound 
states beyond ( ^.7| ), i.e. as indicating that there is a velocity bound (a "speed limit") for the mesons. 



One can also turn (2~7) around and infer that for any large m q and v close to 1, the dissociation 



3 Applying a Lorentz boost to p yields 7 2 (1 + \v 2 )p. Including the (1 + \v 2 ) factor makes this argument reproduce 



the result (2.4), including the weak velocity dependence in the function /, more quantitatively than merely tracking 



the powers of 7. 
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temperature is given by 

T diss = ^(l-*, 2 )3 , (2.8) 

which is consistent with ( p.6|) . Note that the above argument is at best heuristic since Af = 4 SYM 
itself does not contain dynamical quarks and thus genuine mesons do not exist. In the present 
paper, however, we shall see by deriving them from meson dispersion relations that (|2.7|) and (|2.8|) 
are precisely correct in the limit of large quark mass once we introduce fundamentals, and hence 
mesons, into the theory. We shall also find that the more dynamical, albeit heuristic, interpretation 
of ( |2.7D as a velocity beyond which a quark and antiquark do not feel a potential that can bind 
them remains of value. 



3. D3/D7-brane construction of mesons 



In this Section we review the gravity dual description of strongly coupled M = 4 SYM theory with 
gauge group SU(N) coupled toNf<^NAf = 2 hypermultiplets in the fundamental representation 
of SU(N), introduced in gl| and studied in g|, M> M> § 



16 



48, 49, 50, 29, 51, 52 



53. We 
We will 



will first describe the theory at zero temperature and then turn to nonzero temperature 
work in the limit N — > 00, A = gyu-N 00 and Nf finite (in fact Nf = 1). In the deconfined 
strongly coupled plasma that this theory describes, heavy quark mesons exist below a dissociation 

2.166 m q /V\ m,\ 



50, 29 . In Section 5 



temperature that, for mesons at rest, is given by T diss 
we shall calculate the dispersion relations for these mesons, namely the meson spectrum at nonzero 
momentum k and in so doing determine Tdi SS (v) directly, rather than by inference as described in 
Section 2. 



3.1 Zero temperature 

Consider a stack of N coincident D3-branes and Nf coincident D7-branes in 9+1-dimensional 
Minkowski space, which we represent by the array 



D3: 12 3 

D7: 1 2 3 4 5 6 7- 



(3.1) 



which denotes in which of the 9+1 dimensions the D3- and D7-branes are extended, and in which 
they occupy only points. The D3-branes sit at the origin of the 89-plane, with L denoting the 
distance between the D3- and the D7-branes in the 89-plane. Without loss of generality (due to 
rotational symmetry in 89-plane), we can take the D7-branes to be at x§ = L, xg = 0. This is 
a stable configuration and preserves one quarter of the total number of supersymmetries, meaning 
that it describes an M = 2 supersymmetric gauge theory as we now sketch [41 . 

The open string sector of the system contains 3-3 strings, both of whose ends lie on one of the 
N D3-branes, 7-7 strings ending on Nf D7-branes, and 3-7 and 7-3 strings stretching between D3- 
and D7-branes. In the low energy limit 



a 



0. 



L 

2ira' 



finite, 



(3.2) 
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all the stringy modes decouple except for: (i) the lightest modes of the 3-3 strings, which give rise to 
an SU(N) M = 4 SYM theory in 3+1-dimensional Minkowski space; (ii) the lightest modes of the 
3-7 and 7-3 strings, which give rise to Nf hypermultiplets in the TV = 2 gauge theory transforming 
under the fundamental representation of SU(N). The whole theory thus has M = 2 supersymmetry. 
We will call Nf hypermultiplets quarks below even though they contain both fermions and bosons. 
The mass of the quarks is given by 

m « = 2^ ' (3 ' 3) 
where l/(2na') is the tension of the strings. 
In the limit 

N -> oo, Nf = finite, A = gl M N > 1 , (3.4) 



the above gauge theory has a gravity description J41| in terms of D7-branes in the near-horizon 



geometry of the D3-branes, which is AdS$ x S$ with a metric 

ds 2 = (-dt 2 + dx\ + dx\ + dxf) + ^;dr 2 + R 2 dtt 2 5 

r 2 R 2 9 

= (-dt 2 + dx\ + dx\ + dxfj + — dx l > ( 3 - 5 ) 

where r 2 = Y^h=a x \ an d ^^i i s the metric on a 5-sphere. R is the curvature radius of AdS and is 
related to the Yang-Mills theory 't Hooft coupling by 

= VX . (3.6) 

a' 

The string coupling constant g s is related to the gauge theory parameters by 

4tt# s = g 2 YM = , (3.7) 

where g\ M is defined according to standard field theory conventions and is twice as large as the 
Yang-Mills gauge coupling defined according to standard string theory conventions. In this zero 
temperature setting, the embedding of the D7-branes in the AdS 5 x S5 geometry ( |3.5|) can be read 
directly from ( |3.1| ). The D7-brane worldvolumes fill the (t, x{) coordinates, with i — 1, ... ,7, and 
are located at the point x§ = L, xg = in the 89-plane. Since Nf remains finite in the large N 
limit, the gravitational back-reaction of the D7-branes on the spacetime of the D3-branes ( |3.5| ) may 
be neglected. 

The dictionary between the gauge theory and its dual gravity description can thus be sum- 
marized as follows. On the gauge theory side we have two sectors: excitations involving adjoint 
degrees of freedom only and excitations involving the fundamentals. The first type of excitations 
correspond to closed strings in AdS$ x S5 as in the standard AdS/CFT story. The second type is 
described by open strings ending on the D7-branes 4 . In particular, the low-lying (in a sense that we 



We will not consider baryons in this paper. 
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shall define later) meson spectrum of the gauge theory can be described by fluctuations of Xg g and 
gauge fields on D7-branes. We shall focus on the fluctuations of xs,9 on the D7-brane (equivalently, 
the fluctuations of the position of the D7-brane in the (xg,xg) plane) which describe scalar mesons. 
There are also gauge fields localized within the D7-brane, and their fluctuations describe vector 
mesons. The description of the vector mesons is expected to be similar to that of the scalar mesons. 
We shall limit our presentation entirely to the scalar mesons. We shall take Nf = 1, meaning that 
the gauge theory is specified by the parameters N, A and m q which are related to their counterparts 
in the dual gravity theory by (|3.6| ), ( |3.7|) and (|3.3|) . We see that the N — > oo limit corresponds to 
g s — > 0, making the string theory weakly coupled. Considering the theory with the parameter A 
taken to oo corresponds to taking the string tension to infinity. These limits justify the use of the 
classical gravity approximation in which we consider strings moving in a background spacetime. 

For later generalization to finite temperature, it is convenient to describe the D7-brane in a 
coordinate system which makes the symmetries of its embedding more manifest. We split the M 6 
factor in the last term of ( |3.5|) into M 4 x M 2 (i.e. parts longitudinal and transverse to the D7-brane) 
and express them in terms of polar coordinates respectively. More explicitly, 

r 2 = p 2 + y 2 , p 2 = x\ + xl + xl + x 2 , y 2 = x\ + x|, 

x 8 = ycos(f>, x 9 = ysm<p. (3.8) 

The metric (13.51) then becomes 



ds 2 = (-dt 2 + dx 2 ) + -^-5 (dp 2 + p 2 dn 2 3 + dy 2 + y 2 d<j> 2 ) . (3.9) 

The D7-brane now covers (t, x) = (t, xi, X2, £3, p, Q3) and sits at y = L and = 0. Note that in the 
radial direction the D7-brane extends from p = 0, at which the size of the three-sphere ^3 becomes 
zero, to p = 00. The point p = corresponds to r = L. 

We now briefly describe how to find the low-lying meson spectrum described by the fluctuations 
of Xg g. The action of the D7-brane is given by the Dirac-Born-Infeld action 



Sdy = -Pi J Ay-det/iy , (3.10) 

where the (with i = 0, 1, . . . , 7) denote the worldvolume coordinates of the D7 brane and is 
the induced metric in the worldvolume 

The value of the D7-brane tension, p,-j = (27T)^ e g~ 1 a'~ A , will play no role in our considerations. 
The spacetime metric is given by (|3.9| ) and X M (£) describe the embedding of the D7-brane, 
where p, runs through all spacetime coordinates. The action ( |3.10| ) is invariant under the coordinate 
transformations £ — > £'(£)• We can use this freedom to set £ l = (t,x, p,Q 3 ), and the embedding 
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described below equation ( |3.9| ) then corresponds to the following solution to the equations of motion 
of (CT: 

y{£) = L, 0(0 = or x 8 (£) = £, x 9 (£) = . (3.12) 
To find the meson spectrum corresponding to the fluctuations of the brane position, we let 



x 8 = L + 27raVi(0 , x 9 = + 27ra / V>2(£); 
and expand the action ( |3.10| ) to quadratic order in ^ 12 , obtaining 



S D7 ~ fj, 7 I d^p 3 [ -1 - ^(2ira'R) 2 



p 2 + L 2 



(3.13) 



(3.14) 



In (|3.14j) , hij denotes the induced metric on the D7-brane for the embedding Q3.12 ) in the absence 
of any fluctuations, i.e. 



ds 2 = h^icde = p -^- (-dt 2 + dx 2 ) + (dp 2 + P 2 dnl) . 



R 2 



p 2 + L 2 



(3.15) 



Note that when L = 0, the above metric reduces to Ad S§ x S 3 , reflecting the fact that in the 
massless quark limit the Yang-Mills theory is conformally invariant in the large N/Nf limit. 

The equation of motion following from ( |3.14| ) is 



R 4 



1 d 



(p 2 + L 2 ) 2 



p 3 dp \ dp 



d 



d a d a ^ -\- — — I p d — $ + -V^ = , 



(3.16) 



where ip denotes either ipi or ip2, where a = ... 3, and where V 2 denotes the Laplacian operator on 
the unit S* 3 . Eq. ( p,16|) can be solved exactly and normalizable solutions have a discrete spectrum. 
It was found in El that the four dimensional mass spectrum is given by 



m n i 



~7T 



y/(n + l + l)(n + l + 2), 



n = 0,l,... , Z = 0,1,... , 



(3.17) 



with degeneracy (£ + l) 2 , where I is the angular momentum on S* 3 . The (£ + l) 2 degeneracy is 
understood in the field theory as arising because the scalar mesons are in the (£/2,£/2) representa- 
tion of a global SU(2) x SU(2) symmetry corresponding to rotations in the S 3 in the dual gravity 
theory EJ . 



The mass scale appearing in (|3.17|) can also be deduced without calculation via a scaling 
argument. Letting 



R 2 

T ' 



x -> —x, p -> Lp, 



the metric ( |3.15| ) can be solely expressed in terms of dimensionless quantities: 

^ = (p 2 + 1) (-dt 2 + dx 2 ) + - J^ (dp 2 + p 2 dn 2 3 ) . 



(3.18) 



(3.19) 
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Thus, the mass scale for the mesonic fluctuations must be 



Ms w = i^- (3 ' 20) 

as is indeed apparent in the explicit result ( |3.17|) . We see that the mesons are very tightly bound 
in the large A limit with a mass M that is parametrically smaller than the rest mass of a separated 
quark and antiquark, 2m q . This means that the binding energy is ~ —2m q . From this fact and the 
Coulomb potential (|2.1| ), one can also estimate that the size of the bound states is parametrically 



of order - ;/ 



1 



Finally, we can now explain the sense in which our analysis is limited to low-lying mesons. 
We are only analyzing those scalar mesons whose mass is of order M. There are other, stringy, 
excitations in the theory with meson quantum numbers whose masses are of order L/(R\/c/) ~ 
MX ' ~ m n /A 1 / 4 and of order L/a' ~ MX 1 ' 2 ~ m q [42]. They are parametrically heavier than 



the mesons we analyze, and can be neglected in the large A limit even though those with masses 
~ m g /A 1//4 are also tightly bound, since their masses are also parametrically small compared to m q . 
In Section 5, we shall see again in a different way that our analysis of the dispersion relations for 
the mesons with masses ~ m q /\/X that we focus on is controlled by the smallness of 1/A 1//4 . 



3.2 Nonzero Temperature 

We now put the Yang-Mills theory at nonzero temperature, in which case the AdS§ part of the 
metric ( |3.5| ) is replaced by the metric of an AdS Schwarzschild black hole 

r 2 1 

ds 2 = -f(r)dt 2 + —rdx 2 + -rr-Ar 2 + R 2 dQ\ , (3.21) 
R j{r) 

fir) (3-22) 

The temperature of the gauge theory is equal to the Hawking temperature of the black hole, which 
is 

T — —jyj ■ (3.23) 

This is the one addition at nonzero temperature to the dictionary that relates the parameters of 
the (now hot) gauge theory to those of its dual gravity description. 

At nonzero temperature, the embedding of the D7-brane is modified because the D7-brane now 
feels a gravitational attraction due to the presence of the black hole. To find the embedding, it is 
convenient to use coordinates which are analogous to those in (|3.9| ). For this purpose, we introduce 
a new radial coordinate u defined by 



dr 2 R 2 du 2 . o 1 



f(r) u 2 2 



i.e. u 2 = - [r 2 + Wr 4 - r 4 , (3.24) 
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in terms of which ( |3.21 ) can then be written as 

ds 2 = -fdt 2 + ^-dx 2 + ^(du 2 + u 2 dnl) (3.25) 
R 2 u 

= -fdt 2 + ^dx 2 + ^ (dp 2 + p 2 dn 2 3 + dy 2 + y 2 d<p 2 ) . (3.26) 

As in ( p.9|) , we have split the last term of ( [3.25|) in terms of polar coordinates on M 4 x M 2 , with 

u 2 = y 2 + p 2 . (3.27) 
In ( |3.25| ) and Q3.26| ), / and r should now be considered as functions of u, 



4 („A „4/a\2 



,2 



Au 2 u 2 R 2 (u A + 7o/4j 



r 



r 4 /4) 



In terms of w, the horizon is now at Uo = ^7=. 

The D7-brane again covers £ l = (t, x, p, Jl 3 ) and its embedding y(£),<j>(£) in the (y, <fi) plane will 
again be determined by extremizing the Dirac-Born-Infeld action (|3.10|) . Because of the rotational 
symmetry in the <fi direction, we can choose </>(£) = 0. Because of the translational symmetry in the 
(t,x) directions and the rotational symmetry in S 3 , y can depend on p only. Thus, the embedding 
is fully specified by a single function y(p). The induced metric on the D7-brane worldvolume can 
be written in terms of this function as 

hudgd? = -f(u)dt 2 + ^dx 2 + ^ ((1 + y'(p) 2 )dp 2 + p 2 dnj) , (3.29) 

where u in (13.271 ) and hence f(u) are functions of p and y(p). Substituting ( |3.29| ) into ( |3.10| ), one 
finds 

SW X I dp ^pf ( 16 0^) 8 " *) ' /TT¥W2 ' (3 ' 30) 
which leads to the equation of motion 



y" | V | 8r 8 jpy'-y) 
1 + y' 2 p u 2 16u 8 — r 



72 + ^ + ^7^3 = ° ( 3 - 31 ) 



for y(p), where u 2 (p) = p 2 + y 2 (p). 

To solve ( |3.31| ) one imposes the boundary condition that y —>■ L as p — > oo, and that the 
induced metric (|3.29|) is non-singular everywhere. L determines the bare quark mass as in (|3.3j ). It 
is convenient to introduce a parameter 

2 2 \ T 2 2 T 2 

00 " L 2 2L 2 8m 2 2M 2 ' 1 } 
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Figure 1: Some possible D7-brane embeddings y(p). The quark mass to temperature ratio is determined 
by y(oo) = L. Specifically, V8m q /(T\/ r \) = y(oo)/uo = 1/y/Coo- The top three curves are Minkowski 
embeddings, with y{p) extending from p = to p = oo. The bottom three curves are black hole embeddings, 
in which the D7-brane begins at the black hole horizon at y 2 + p 2 = Uq. The middle curve is the critical 
embedding. The seven curves, ordered from top to bottom as they occur in the left part of the figure, 
are drawn for temperatures specified by = 0.249, 0.471, 0.5865, 0.5948, 0.5863, 0.647 and 1.656. Note 
that the = 0.5863 black hole embedding crosses both the eoo = 0.5948 critical embedding and the 
£oo = 0.5865 Minkowski embedding. 



where we have used (|3.23| ) and ( [3.20]) . Because J\ — 4 SYM is scale invariant before introducing the 



massive fundamentals, meaning that all dimensionful quantities must be proportional to appropriate 
powers of T, when we introduce the fundamentals the only way in which the quark mass m q can 
enter is through the dimensionless ratio m q /T. Scale invariance alone does not require that this 
ratio be accompanied by a y/\, but it is easy to see that, after rescaling to dimensionless variables 
as in (|3.18|) , the only combination of parameters that enters (|3.31|) is e^. The small e M regime can 
equally well be thought of as a low temperature regime or a heavy quark regime. In the remainder 
of this section, we shall imagine m q as fixed and describe the physics as a function of varying T, 
i.e. varying horizon radius r . 

The equation of motion ( [3.3 1|) that specifies the D7-brane embedding can be solved numerically. 



Upon so doing, one finds that there are three types of solutions with different topology [ID], ^DJ : 

• Minkowski embeddings: The D7-brane extends all the way to p = with y(0) > uq = ^= 
(see e.g. the upper three curves in Fig. |l|). In order for the solution to be regular one needs 
y'{0) = 0. This gives rise to a one-parameter family of solutions parameterized by y(0). The 
topology of the brane is M 1 ' 7 . 
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Critical embedding: The D7-brane just touches the horizon, i.e. y(0) = u (see e.g. the middle 
curve in Fig. [TJ). The worldvolume metric is singular at the point where the D7-brane touches 
the horizon. 

Black hole embeddings: The D7-brane ends on the horizon w = r /\/2 at some p > (see 
e.g. the lower three curves in Fig. [I]). The topology of the D7-brane is then M 1,4 x S 3 . 

It turns out |K| ^9| that Minkowski embeddings that begin at p = with y close to r /v^2, 



just above the critical embedding, can cross the critical embedding, ending up at p — >• oo with 
y(oo) just below that for the critical embedding. Similarly, embeddings that begin just below the 
critical embedding can end up just above it. Furthermore, those embeddings that begin even closer 
to the critical embedding can cross it more than once. This means that there is a range of values 
around the critical = 0.5948 for which there are three or more embeddings for each value of 
eoo. At low temperatures (precisely, for < 0.5834) this does not occur: there is only a single 
Minkowski embedding for each value of e ro . At high temperatures (precisely, for 6^ > 0.5955) there 
is only a single black hole embedding per value of eoo. In the intermediate range of temperatures 
0.5834 < < 0.5955, one needs to compare the free energy of each of the three or more different 
D7-brane embeddings that have the same value of Coo to determine which is favored. One finds 
that there is a first order phase transition at a temperature T c at which eoo = 0.5863, where the 
favored embedding jumps discontinuously from a Minkowski embedding to a black hole embedding 
as a function of increasing temperature [[50], p9]. 5 



As we shall study in detail in Section 4, fluctuations about a Minkowski embedding describe a 
discrete meson spectrum with a mass gap of order O(M). In contrast, fluctuations about a black 
hole embedding yield a continuous spectrum |5l], |29|. A natural interpretation of the first order 
transition is that T c = T diss , the temperature above which the mesons dissociate |)0j p9| . It is 
interesting, and quite unlike what is expected in QCD, that all the mesons described by the zero 



temperature spectrum (|3.17| ) dissociate at the same temperature. This is presumably related to 
the fact that the mesons are so tightly bound, again unlike in QCD. We shall therefore focus on 
the velocity-dependence of the meson spectrum at nonzero temperature — in other words, the 
meson dispersion relations first studied in [p9| . As we have explained in Section 1, the velocity- 



dependence is currently inaccessible to lattice QCD calculations. Hence, even qualitative results are 
sorely needed. Furthermore, inferences drawn from a previous calculation of the potential between 
a moving quark-antiquark pair lead to a velocity-scaling (j2.6|) of Tdi SS that has a simple physical 



interpretation, which suggests that it could be applicable in varied theories p2| . We shall see this 
velocity dependence emerge from the meson dispersion relations in Section 5. 

It is interesting to return to the qualitative estimate of T^ iss obtained from the static quark- 
antiquark potential in Section 2, and see how it compares to the Tdi SS = T c obtained from the analysis 
of the mesons themselves. Equating the size of a meson with binding energy 2m q , determined by 
the zero-temperature potential ( ^.1| ), with the screening length L s = L c = 0.24/T, determined by 



5 The critical embedding occurs at an = 0.5948 which is greater than the e M at which the first order phase 
transition occurs, meaning that at e w = 0.5948 there is a black hole embedding that has a lower free energy than 
the critical embedding. 
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Figure 2: (determined by the embedding y at infinity) versus e (determined either by y(0), for 
Minkowski embeddings with e < 1, or by where the embedding intersects the horizon, for e > 1). The 
right panel zooms in on the vicinity of the critical embedding at e = 1. The stable embeddings and the 
first order phase transition are indicated by the thick curve; the metastable embeddings are indicated by 
the thin curves. 



the potential ( |2.2| ) at nonzero temperature, yields the estimate that Tdi SS should be ~ 2.1m q / V\. 
This is in surprisingly good agreement with T diss = y/8e^m q / V\ = 2.166 m q /VX for Coo = 0.5863. 

In subsequent sections, we shall derive the dispersion relations for mesons at T < T^ [ss . We 
close this section by introducing some new notation that simplifies the analysis of the Minkowski 
embedding of the D7-brane, whose fluctuations we shall be treating. We first introduce parameters 



2 2 

u 2 _ r 2 



V(U) . £ = 75 = 7T72 ■ ( 3 - 33 ) 



T 2 9T 2 



For Minkowski embeddings, e takes value in the range [0,1], with e = corresponding to zero 
temperature, and e = 1 to the critical embedding. Although that we introduced earlier has the 
advantage of being directly related to the fundamental parameters of the theory according to ( |3.32| ), 



the new parameter has the advantage that there is only one embedding for each value of e. And, 
e will turn out to be convenient for analyzing the equations of motion ( |3.3 1| ) and the fluctuations 
on D7-branes. When too <C 1, e ~ t^. A full analytic relation between e and is not known, 
but given an e one can readily find the corresponding numerically. For example, at T = T c , 
e = 0.756 and e ro = 0.586 while for the critical embedding, e = 1 and e M = = 0.5948. We depict 
the relation between and e in Fig. |2|. In order to make this figure, for the black hole embeddings 
we have defined £ — 1/ sin 2 6 where 6 is the angle in the (y, p) plane of Fig. [I] at the point at which 
the black hole embedding y(p) intersects the black hole horizon y 2 + p 2 = Uq. That is, 1 < e < oo 
parametrizes black hole embeddings which begin at different points along the black hole horizon. 
The seven embeddings in Fig. [I] have e = 0.25, 0.5, 0.756, 1.00, 1.13, 1.41 and 4.35, from top to 
bottom as they are ordered on the left, i.e. at the tip of the D7-brane at y = for the Minkowski 
embeddings and at the horizon for the black hole embeddings. 
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Finally, it will also prove convenient to introduce dimensionless coordinates by a rescaling 
according to 

R 2 R 2 

t — > —t, Xi — > —Xi, p — > L p, y — > L y, (3.34) 

after which the spacetime metric becomes 

ds 2 1 

— r = G^dx^dx" = -f(u)dt 2 + riufdx 2 + (dp 2 + p 2 dQ\ + dy 2 + y 2 dcf> 2 ) (3.35) 
R z u 

and the induced metric becomes 

r7s 2 1 

= hijdCdg = -f(u)dt 2 + r 2 dx 2 + — ((1 + y\pf)dp 2 + p 2 dQj) (3.36) 
R z u z 

with 

u 2 = y 2 + p 2 , f(u)= r 2 (u)=u 2 + £ - 2 , (3.37) 

u z \ur + e A ) u 

where both G^ u and are now dimensionless. The equation of motion for y(p) becomes 



1 + y' 2 p u 2 \e~ 4 u 8 — 1 
with the boundary conditions 

2/(0) = 1, ?/(0) = . (3.39) 

This form of the equations of motion that determine the embedding y(p) will be useful in subsequent 
sections. 



4. Meson Fluctuations at Nonzero Temperature 

In this section we derive linearized equations of motion that describe the small fluctuations of 
the D7-brane position. A version of these equations have been derived and solved numerically by 



various authors (see e.g. [43], |4], |5(J, ^9|). Here we will rederive the equations in a different form 
by choosing the worldvolume fields parameterizing the fluctuations in a more geometric way. The 
new approach gives a nice geometric interpretation for the embedding and small fluctuations. It 
also simplifies the equations dramatically, which will enable us to extract analytic information for 
the meson dispersion relations in the next section. We present the main ideas and results in this 
Section but we leave technical details to Appendix A. In that Appendix, we shall also present a 
general discussion of the fluctuations of a brane embedded in any curved spacetime. 

The action for small perturbations of the D7-brane location can be obtained by inserting 

x"(0 = XS(C) + sx^c) (4.i) 
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into the D-brane action ( 3.10 ) and ( 3.11 ), where Xtf (£) denotes the background solution that de- 



scribes the embedding in the absence of fluctuations, and SX^ describes small fluctuations transverse 



to the brane. For the D7-brane under consideration, in the coordinates used in Q3.35| ) the general 
expression (|4.1|) becomes 

y(Z) = yo(p) + 6y(0, 0(0 = 50(0 (4.2) 

with yo(p) the embedding solution obtained by solving (|3.38|) . The choice of the worldvolume fields 
5y, 5<p is clearly far from unique. Any two independent functions of 5y and 50 will also do. (This 
freedom corresponds to the freedom to choose different coordinates for the 10-dimensional space 
within which the D7-brane is embedded.) In fact, it is awkward to use by and 50 as worldvolume 
fields since they are differences in coordinates and thus do not transform nicely under coordinate 
changes. Using them obscures the geometric interpretation of the equations. Below we will adopt 
a coordinate system which makes the geometric interpretation manifest. Since our discussion is 
rather general, not specific to the particular system under consideration, we will describe it initially 
using general language. 

Consider a point Ao(0 on the brane. The tangent space at Xq perpendicular to the D7-brane 
is a two-dimensional subspace Vq spanned by unit vectors nl, which are orthogonal to the branes, 
i.e. 

<^(^-T-y'o(p)(i-Y (4-3) 



dy J \dp 
^ocll-V. (4.4) 



Any vector rj' 1 in Vq can be written as 

rf = Xl n^ + x 2 ^2 • ( 4 - 5 ) 

We can then establish a map from {xiiXz) to small perturbations SX^ in (|4.1| ) by shooting out 
geodesies of unit affine parameter from X with tangent r] fl . Such a map should be one-to-one for 
Xii X2 sufficiently small. Clearly \i and X2 behave like scalars under coordinate changes and we will 
use them as the worldvolume fields parameterizing small fluctuations of the position of the brane. 
By solving the geodesic equation, 5X^ can be expressed in terms of xi,2 as 

sx* = 1 f- V + • • • , (4.6) 

where F^p are the Christoffel symbols of the 10-dimensional metric. Note that the choice of xi,2 
is not unique. There is in fact an 5*0(2) "gauge" symmetry under which xi,2 transform as a 
vector, since one can make different choices of basis vectors n\, n 2 that are related by a local SO (2) 
transformation. 

We now insert ( [4.6[ ) and ( |4.1| ) into the Dirac-Born-Infeld action ( |3.10| ) and, after some algebra 
discussed further in Appendix A, we find that the equations of motion satisfied by Xq (i.e. which 
determine the embedding in the absence of fluctuations) can be written as 

K s = 0, (4.7) 
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and the quadratic action for small fluctuations Xi,2 about X takes the form 



Sdj = ViR* j (PtyZ-dethij (-^DiXsD'xs - \xsXt (~K sij Kf + i? si ^) J 
where s, t — 1,2 and where we have defined the following quantities: 



(4.8) 



h 



'j 



K 

DiXs 



G^diXQdjXQ , 
= diXs + U ist xt , 



R 



sijt 



•tfr&diXRdiXZFL 



K W 

v Sl'J 1 b 1 



U, 



ist 



n 



(4.9) 
(4.10) 
(4.11) 



Note that hij is the induced metric on the brane and i,j are raised by h l K R af iuf3 is the Riemann 
tensor for the 10-dimensional spacetime. K S ij is the extrinsic curvature of the brane along the 
direction n^. Ui S t (which is antisymmetric in s,t) is an SO (2) connection for the SO (2) gauge 
symmetry and D$ is the corresponding covariant derivative. We see that the embedding equations 
of motion ( |4.7| ) have a very simple geometric interpretation as requiring that the trace of the 
extrinsic curvature in each orthogonal direction has to vanish, which is what we expect since this 
is equivalent to the statement that the volume of the D7-brane is extremal. 

The symmetries of the D7-brane embedding that we are analyzing allow us to further simplify 
the action (4.8). Because in (|4.4j) is proportional to a Killing vector and is hypersurface or- 



thogonal, Z7ji2 and K 2 ij vanish identically. (See Appendix A for a proof, and for the definition of 
hypersurface orthogonal.) With K 2 = satisfied as an identity, the remaining equation of motion 
specifying the embedding, namely K\ = 0, is then equivalent to the equation of motion for y that 
we derived in Section 3, namely Eq. (|3.38|) . After some further algebra (see Appendix A) we find 
that the action (|4.8|) for small fluctuations reduces to 



S D7 = fi 7 R 8 J d^-dethij ^~(dxi 



;(,dx2 
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where we have defined 



-R2112 

R11 
R22 
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(4.14) 



and where R is the Ricci scalar for the 10-dimensional spacetime while ( 8 'R is the Ricci scalar for the 
induced metric on the D7 brane. The background metric is given by ( p. 35 ). The "masses" 
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Figure 3: The squared "masses" of the two orthonormal geometric modes of the D7-brane fluctuations 
for Minkowski embeddings (left panel) and black hole embeddings (right panel). In each figure, ml (m^) is 
plotted as a solid (dashed) line for three values of e^. The Minkowski embeddings have eoo = 0.587, 0.471 
and 0.249 (top to bottom) and the black hole embeddings have = 1.656, 0.647 and 0.586 (again top to 
bottom, this time with temperature increasing from top to bottom.) The Minkowski embedding is plotted 
as a function of p and the black hole embedding as a function of u with the horizon on the left at u = 1. 



m\ and m\ are nontrivial functions of p. Since the worldvolume metric is regular for Minkowski 
embeddings, they are well denned for p e [0, oo). 

Our result in the form (|4.8|) is very general, applicable to the embedding of any codimension- 
two branes in any spacetime geometry. For example, we can apply it to the embedding of D7-branes 
at zero temperature given by ( |3.12| ) and learn that the meson fluctuations at zero temperature are 
described by ( [4.12j ) with 

3p 2 + 4 



m-, 



ml 



l + p 2 



(4.15) 



and with hij in ( |4.12|) given by (p.15 ). It is also straightforward to check that equations of motion 
derived from ( |4.12[ ) with ( [4.15 ) and given by ( 3.15 ) are equivalent to ( 3.16| ). At zero tempera- 
ture, ( 3.14 ) and ( |3.15|) are already simple enough and the formalism we have described here does 
not gain us further advantage. However, at nonzero temperature the equations of motion obtained 
from (|4.12|) yield both technical and conceptual simplification. In Section 5 we shall use the for- 
malism that we have developed to obtain the dispersion relations at large momentum analytically. 

Before turning to the dispersion relations, we plot the "masses" ml and m\ for various D7- 
brane embeddings at nonzero temperature in Fig. |3|. Using a numerical solution for y(p), it is 
straightforward to evaluate (|4.13| ), obtaining the masses in the figure. For the black hole embeddings, 
the D7-brane begins at the black hole horizon at u = 1 rather than at p = 0, see Fig. p], making it 
more convenient to plot the masses as a function of u rather than p. We can infer several important 
features from the masses plotted in Fig. |3|. As p — > oo, both ml and m\ approach —3 for all the 
embeddings. This implies that xi an d Xi couple to boundary operators of dimension 3, as shown 
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in [51] by explicit construction of the operators in the boundary theory which map onto X i an d X2- 
As e — > 1 from below for the Minkowski embeddings (from above for the black hole embeddings), 
the behavior of m 2 at the tip of the D7-brane at p = (at u — 1) becomes singular, diverging to 
minus infinity This is a reflection of the curvature divergence at the tip of the critical embedding 
at p = (u = 1). 

We have referred to mf and m 2 as "masses" in quotes because the equations of motion obtained 



by straightforward variation of the action (|4.12 ) in which they arise yields 




hh t >d jX .)-m 2 a x. = 0, s = l,2 (4.16) 



with h = dethij, which is a Klein-Gordon equation in a curved spacetime with spatially varying 
"masses". If we could cast the equations of motion in such a way that they take the form of a 
Schrodinger equation with some potential, this would make it possible to infer qualitative implica- 
tions for the nature of the meson spectrum immediately via physical intuition, which is not possible 
to do by inspection of the curves in Fig. |3|. To achieve this, we recast the equations of motion as 
follows. We introduce a "tortoise coordinate" z defined by 

^ 2 = ^(1+?/o(p)W, (4-17) 
in terms of which the induced metric on the brane takes the simple form 

^f|l = f(-dt 2 + dz 2 ) + r 2 (u)dx 2 + £dn 2 3 . (4.18) 
R 2 u 

(We choose the additive constant in the definition of z so that z = at p = 0.) Then, we seek 
solutions to the equations of motion ( 4.16|) that separate according to the ansatz 



Xs = ^ e~^ t+rk - £ Y emrh (Q 3 ) (4.19) 

with 



i 

2 , 3 

'rp\i 



f j w ■ (4 - 20) 

Such a solution is the wave function for a scalar meson of type s = 1 or s = 2 with energy uj and 
wave vector k (note the plane wave form for the dependence on (3+l)-dimensional Minkowski space 
coordinates) and with quantum numbers £, m and m specifying the angular momentum spherical 
harmonic on the "internal" three-sphere. (Rotation symmetry of the three-sphere guarantees that 
the quantum numbers m and fh will not appear in any equations.) The ip s (z) that we must solve 
for are the wave functions of the meson states in the fifth dimension. 

The reasons for the introduction of the tortoise coordinate z and the ansatz ( 4.19Q for the form 



of the solution become apparent when we discover that the equations of motion ( 4.16Q now take the 
Schrodinger form 

d 2 



dz- 



i) s + V s (k, z)il) s = uj 2 ^ s , (4.21) 
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Figure 4: Potentials V s {z) for Minkowski embeddings at various temperatures, all with k = £ = 0. The 
left (right) panel is for s = 1 (s = 2). In each panel, the potentials are drawn for eoo = 0.249, 0.471, 0.586 
and 0.5948, with the potential widening as the critical embedding is approached, i.e. as is increased. 
The Eqo = 0.586 potential is that for the Minkowski embedding at the first order transition; the widest 
potential shown describes the fluctuations of a metastable Minkowski embedding very close to the critical 
embedding. The potential becomes infinitely wide as the critical embedding is approached, but it does so 
only logarithmically in 6^ — 6qo • Note that the tip of the D7-brane is at z — 0, on the left side of the figure, 



whereas p = oo has been mapped to a finite value of the tortoise coordinate z 
the "wall" on the right side of each of the potentials in the figure. 



corresponding to 



V s (k, z) = fL r + fm l + ^-+ ^UfL . (4.22) 



with potentials for each value of k — \k\ and for each of the two scalar mesons labelled by s = 1, 2 
given by 

Z r z p z 

Here, the prime denotes differentiation with respect to z. Recall that u 2 = p 2 + Uq(p) and it should 
be understood that p, u, and yo are all functions of the tortoise coordinate z. In Figs. f| and |, we 
provide plots of V s (z) with k = I = for s — 1,2 and for Minkowski (Fig. f|)) and black hole (Fig. ||) 
embeddings at various temperatures. With the tortoise coordinate z denned as we have described, 
in a Minkowski embedding z extends from z — 0, which corresponds to the tip of the D7-brane, to 



which corresponds to p = oo. Here, u(p) and f(u) are given in ( |3.37p . This defines the width of the 



potentials for the Minkowski embeddings shown in Fig. [|, which get wider and wider as the critical 
embedding is approached. 

If we used the same tortoise coordinate for the black hole embeddings, the lower limit of the 
integral (|4.23|) is then the p at which y(p) intersects the horizon and f(u) vanishes, making the 
integral divergent. This means that p = oo is mapped onto z = oo for black hole embeddings. 
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Figure 5: Potentials V s {zbh) for black hole embeddings at various temperatures, all with k = I = 0. The 
left (right) panel is for s = 1 (s = 2). In each panel, the potentials are drawn for eoo = 3584., 0.647, 0.586, 
0.586, 0.5940 and 0.5948, from narrower to wider, with the potential widening as the critical embedding 
is approached from the right along the curve in Fig. ^. Note that Zbh is defined such that the horizon is at 
Z\,h = oo, and p = oo is at Zbh = 0. The narrower (wider) of the two potentials with 6^ = 0.586 is that for 
the stable (unstable) black hole embedding: at this e^, there is a first order transition (see Fig. §) between 
the stable Minkowski embedding (whose potential is found in Fig. and the stable black hole embedding. 
The potentials at = 0.5940 and 0.5948 describe fluctuations of metastable black hole embeddings, with 
the latter being a black hole embedding very close to the critical embedding. 



It is more convenient to define Zbh by first choosing the integration constant such that p = oo 
corresponds to Zbh — 0, and then multiplying by -1. This is the tortoise coordinate that we have 
used in Fig. || 

The qualitative implications for the meson spectrum can be inferred immediately from Figs. ^ 
and [|, since we have intuition for solutions of the Schrodinger equation. We can see immediately that 
the Minkowski embeddings all have a discrete spectrum of meson excitations, while the fluctuations 
of the black hole embeddings all have continuous spectra. This justifies the identification of the 
first order phase transition from Minkowski to black hole embeddings that we described in Section 
3 as the transition at which mesons dissociate. 



Other phenomena that are discussed quantitatively in [29, 50, 51 1 can be inferred qualitatively 



directly from the potentials in Figs. |] and |5|. For example we see from the left panel in Fig. that, in 
addition to the continuous spectrum characteristic of all black hole embeddings, those embeddings 
that are close to the critical embedding will have discrete bound states for the ip\ fluctuations. These 
bound states will always have negative mass-squared, representing an instability. This instability 
arises only in a regime of temperatures at which the black hole embeddings already have a higher 
free energy than the stable Minkowski embedding, that is, at temperatures below the first order 



transition p9| . They therefore represent an instability of the branch of the spectrum that was 
already metastable. Similarly, the left panel of Fig. f| shows that Minkowski embeddings close to 
the critical embedding also have negative mass-squared bound states; again, this instability only 
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occurs for embeddings that were already only metastable [29]. We see from the right panel of Fig. |5] 
that resonances may also occur in the ip2 channel for the black hole embedding. They are interpreted 
as quasi-normal modes; close to the transition these resonances become more well defined and may 



be interpreted as quasi-particle meson excitations pi |T 



5. Dispersion relations 



We have now laid the groundwork needed to evaluate the dispersion relations for the ipi and ip2 
scalar mesons, corresponding in the gravity dual to fluctuations of the position of the D7 brane. 
These fluctuations are governed by ( 4.21|) , which are Schrodinger equations with the potentials 



Vx(k,z) or V 2 (k,z) given by (|4.22| ) and ( |4.13| ) and depicted in Fig. £|. The eigenvalues of these 
Schrodinger equations are u 2 for the mesons. So, it is now a straightforward numerical task to 
find the square root of the eigenvalues of the Schrodinger equation with, say, potential Vi(k, z), at 
a sequence of values of k. At k = 0, this will reproduce the results that we reviewed in Section 
3.2. As we increase k, we map out the dispersion relation u of each of the ip\ mesons. In Fig. |8| in 
Section |5~3] below, we show the dispersion relations for the ground state ipi meson at several values 



of the temperature. Such dispersion relations have also been obtained numerically in [22]. In order 



to more fully understand the dispersion equations, and their implications, we shall focus first on 
analytic results. The potentials are complicated enough that we do not have analytic solutions for 
the general case. We shall show, however, that in the low temperature and/or the large-fc limit, the 
equations simplify sufficiently that we can find the dispersion relations analytically. It is the large-fc 
limit that is of interest to us, but it is very helpful to begin first at low temperatures, before then 
analyzing the dispersion relations in the large-fc limit at any temperature below the dissociation 
temperature. 

Readers who are only interested in the final results can proceed directly to Section |5.4| , where 
we summarize and discuss our central results for the dispersion relations. 



5.1 Low temperature 

At low temperature, e C 1, the D7-branes are far from the horizon of the black hole. In this regime, 
we can expand various quantities that occur in the potentials ( |4.22| ) as power series in e 2 . We shall 
then be able to determine the dispersion relations analytically to order e 2 in two limits: (i) e — > 
at fixed k, meaning in particular that ek — > 0; and (ii) e — > at fixed, large, ek, meaning that 
k — > oo as e — > 0. 

We begin by seeing how the equation (|3.38|) that determines the embedding y(p) in the absence 
of fluctuations simplifies at small e. Expanding y(p) as a power series in e, one immediately finds 
that y(p) is modified only at order e A , i.e. 

y(p) = 1 + 0{e A ) , (5.1) 

which in turn implies that 

eoo = e (1 + 0(e 4 )) . (5.2) 
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Thus, if we work only to order e 2 , we can treat the embedding as being y(p) = 1, as at zero 
temperature, and can neglect the difference between e and (which is to say the difference between 
y(0) and y(oo)). From ( |3.37|) , then, 

u i = 1 + p i + (e 4 ), f( u )^ u ^l-^ + 0(e^ . (5.3) 

By expanding the curvature invariants in ( |4 . 1 3| ) to order e 2 , we find that 

m? = ml = -i±^ + 0(e*) , (5.4) 
1 + (r 

meaning that to order e 2 the mass terms occurring in ( |3.38| ) are as in (ft. 15 ) at zero temperature. 
Next, we expand the tortoise coordinate Q4.17Q , finding 

z = tan- 1 p + e 2 g{p) + 0{e A ), with g(p) = jL ( 3 tan" 1 p + P / 5 + 3/ f ^ \ . (5.5) 



We can then invert ( |5.5| ) to obtain p in terms of z: 



p = tan^-g 2 ^ V - ; + ... . (5.6) 

cos 2 z 

Using these equations, we find that the potential ( f4.22| ) is given to order 0(e 2 ) by 

V(z) = k 2 + V°(z) - 4e 2 k 2 cos 4 z + e 2 h(z) + 0(e\ e*k 2 ), (5.7) 

where 

V°( z ) = -^--l, with a e = ^ + £(£ + 2) (5.8) 
sin 2z 4 

is the potential at zero temperature, and 

2 sin [2z) 4 
We shall not use the explicit form of h(z) in the following. 

5.1.1 Low temperature at fixed k 

At zero temperature (e = 0), solving the Schrodinger equation (|4.21|) with potential V°(z) yields 
the eigenvalues (and hence the dispersion relations) 

iu 2 -k 2 = m 2 n£ , n=l,2,..., 1 = 0,1,..., (5.10) 

with m n £ given by ( |3.17| ) (after restoring its dimensions). If we work in the limit e — > with k 
fixed, then both the 0(e 2 ) and the 0(e 2 k 2 ) terms that describe the effects of nonzero but small 
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temperature in the potential (|5.8| ) can be treated using quantum mechanical perturbation theory. 
To first order in e 2 , the dispersion relation becomes 

u 2 = v 2 nl k 2 + m 2 nl + e 2 b n£ + 0(e 4 ) (5.11) 

with 

v 2 e = 1- a nt > e 2 , 

a n £ = 4(n, £\ cos 4 z\n, I) , 

b nl = (n,£\h(z)\n,l) , (5.12) 

where \n,£) are the eigenfunctions of the Hamiltonian with the unperturbed potential V° of (|5.8|) , 
with wave functions 



Using the recursion relations for the generalized Gegenbauer polynomials Cn fl54"| , a n £ can be 
evaluated analytically, yielding 

_o (n + 2Z + l)(n + 2Z + 2) (n + l)(n + 2) 

" £ 4(n + / + l/2)(n + / + 3/2) 4(n + Z + 3/2) (n + Z + 5/2) ' 1 ' ' 

So, for the ground state with n = £ = 0, a 00 = 18/15. b n i can be computed numerically, but we 
will not do so here. The dispersion relation fl5.11|) is valid for e 2 < 1 and e 2 k 2 1, meaning that 
at small e it is valid for k 1/e. No matter how small e is, the perturbation theory breaks down 
for k ~ - and ( p.llj ) does not apply. In other words, the low temperature e —> limit and the high 
meson momentum k — >• oo limits do not commute. Even though fl5.ll ) cannot be used to determine 
the meson velocity at large k, it is suggestive. We shall see below that in the large-/? limit, the 
meson velocity is indeed 1 — 0(e 2 ), but the coefficient of e 2 is not given by ( |5.14j ). 

5.1.2 Low temperature at fixed, large, ek 

To explore the behavior of the dispersion relations in the large-/? limit, we now consider the following 
scaling limit 

e -»• 0, k -> oo, with A 2 = k 2 e 2 = finite. (5.15) 



In this limit, the potential (|4.22|) again greatly simplifies and, consistent with (|5.7|) , becomes 

V(z) = k 2 -\ — 1 - 4A 2 cos 4 z . (5.16) 

sin 2 z 

This potential is valid in the limit (|5.15 ) for any value of A, small or large. If A is small, the 
dispersion relation can be determined using perturbation theory as before, yielding ( |5.11| ) without 
the e 2 b n £ term. In order to analyze the large- k regime, we now consider A>1, and seek to evaluate 



- 28 - 



V(z) -k : 




Figure 6: The potential ( |5.16 ) with e = 0.756 and k = 5, 20 and 100. We see that as A = e 2 k 2 increases, 
the minimum of the potential moves towards z = 0, the potential deepens, and the curvature around the 
minimum increases. 



the dispersion relation as an expansion in 1/A. For this purpose, we notice that as A — > oo the 
potential (|5.16|) develops a minimum at 



®£ \ Z 



as depicted in Fig. ||. The curvature about the minimum is V"(zq) oc A 2 . Thus, if we imagine 
watching how the wave function changes as we take the large- A limit, we will see the wave function 
getting more and more tightly localized around the point z = z$ which gets closer and closer to 
z — 0. That is, the wave function will be localized around the tip of the brane z — 0. This motivates 
us to expand the potential around z — 0, getting 

V{z)-h? + l = a l (\ + ^ + ^- + -AA 2 (l-2z 2 + ^ + ..) . (5.18) 



^z 2 ' 3 ' 15 ' "7 ~" \ ' ' 3 
If we now introduce a new variable £ = A^z, the Schrodinger equation ( POT ) becomes 



where 



- d l + p + ) * + Vrp = Erf, (5.19) 

1 



Q 2 = 32, E = j(uj 2 -k 2 + 4A 2 ) , (5.20) 



and V contains only terms that are higher order in 1/A: 

1 

A 
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Thus to leading order in the large A limit, we can drop the V term in (|5.19 ). Upon so doing, and 



using the expression ( |5.£| ) for ag, the equation ( |5.19| ) becomes that of a harmonic oscillator in 4 
dimensions with mass | and frequency fl. This quantum mechanics problem can be solved exactly, 
with wave functions given by 

^ = C 3/2+ ^ +1) (~ f) e^ 2 , (5.22) 

up to a normalization constant, and with eigenvalues given by 

E n = tt(n + 2), n = 0,l,... (5.23) 
In (|5.22 ), Lv is the generalized Laguerre polynomial of order 

u= ^ . (5.24) 

The allowed values of I are determined by the requirement that v must be a non-negative integer. 
The degeneracy of n-th energy level is • Higher order corrections in 1/A can then be obtained 
using perturbation theory. For example, with the next order correction included, the degeneracy 
among states with different i and the same n is lifted and the eigenvalues are given by 

E nl = Q(n + 2) + °-f + 0(1/A 2 ) (5.25) 

with 

Cne = --^n + 2) 2 + -^{i + 2) . (5.26) 



Thus, in the small-e: limit with A fixed and large, we find using ( 5.20|) that the dispersion relation 
becomes 

u 2 ni = (1 - Ae 2 )k 2 + 4V2(n + 2)ke + c ni + 0(l/k) . (5.27) 

Notice that c n e is negative for the ground state, and indeed for any n at sufficiently small I. We learn 
from this calculation that in the large-fc limit, at low temperatures mesons move with a velocity 
given to order e 2 by v = yl — As 2 = 1 — 2e 2 . Recalling that to the order we are working = e, 
this result can be expressed in terms of T, m q and A using (|3.32|) . In the next subsection, we shall 
obtain the meson velocity at large k for any e . 



5.2 Large-A; dispersion relation at generic temperature 

The technique of the previous subsection can be generalized to analyze the dispersion relation in the 
large- k limit at a generic temperature below the dissociation temperature. For general e < 1, one 
again observes that the potential has a sharper and sharper minimum near the tip of the brane z = 
as k becomes larger and larger. Thus, in the large k limit, we only need to solve the Schrodinger 
equation near z — 0. 
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To find the potential V(z) as a power series in z near z = 0, we need to know the solution y(p) 
of ( |3.38D near the tip of the brane at p = 0: 



e 4 (5 + 5e 4 - 3e 8 



_P , yu-ruo -uo , 4 4 

4 -l 3(e 4 -l) 3 ^ ^ ; 



(5.28) 



At small p, using the expansion of y in ( 5.28| ), we find the tortoise coordinate z has the expansion 

(5.29) 



vTT 



l-£ 5 



-p + cv 



Using ( 5.28|) and (|5.29|) in ( |4.22j ), after some algebra we find 



V s {z) = k 2 v 2 + -il 2 e 2 z 2 + (3^ + ...)+^ + lse + 0(z 



(5.30) 



where 



1 



v 

n 2 



l + e 2 ' 

32(l-£ 2 ) 2 (l+£ 4 ) 



-n 2 e 2 



l + e 2 ) 5 

5 - 36e 2 + 28e 4 - 36e 6 + 5e 8 



24(1 + e 2 ) 3 
+ 2) (| + Ae 2 + §e 4 + As 5 + f £ 8 ) - 56e 4 



l + e 



2A3 



72^ = 7w + 



80e 4 



l + e 



2\3 



(5.31) 
(5.32) 

(5.33) 

(5.34) 

(5.35) 



and where oli is given by (|5.8f) . We can understand why the leading difference between the potentials 
V\ and V2 for the mesons ipi and ^2 arises in this approximation in the constant terms 7^ and 72^ 
as follows. We see from (|4 . 2 2|) that the difference between V\ and V2 comes only from m 2 and m|, 
which do not enter multiplied by A; 2 and so cannot affect v , Q 2 or (3^. Furthermore, m 2 and m 2 , 
are curvature invariants, see ( |4.13| ), and must therefore be smooth as p —* because for Minkowski 
embeddings the D7 brane is smooth at p = 0. This means that m\ and m 2 cannot affect the 
coefficient of 1/z 2 in (|5.30| ) . 

We can now obtain the dispersion relations from the Schrodinger equations with potentials 
( |5.30| ) as we did in the previous subsection. After making the rescaling z = /c _1//2 £, the Schrodinger 
equation ( |4.21| ) takes exactly the form ( |5.19| ), with 



E 



1/ 2 
k {uJ 



v 2 k 2 ) 



(5.36) 



31 



where Q and Vq are given by ( p. 31 ) and ( |5.32 ) respectively, and where V s (z) contains only terms 
that are subleading in the 1/k expansion, and is given by 



Thus, we find the large-A; dispersion relation 

2 - k 2 v\ + kQe{n + 2) + d sn£ + C{l/k) 



(5.37) 



(5.38) 



with 



1 



:i+e 2 ) 3 

'5 



~£(£ + 2) (l + 3e 2 + e 4 + 3e 6 + e 8 ) 



4 



- 9e z + 7e q 



9e b + -e l 
4 



(n + 2) 



5Qe 4 



and 



dint + 



80e 4 



1 + e 



2^3 



(5.39) 



(5.40) 



Restoring dimensionful quantities in the dispersion relation (|5.38|) , i.e. undoing (|3.34|) , means 



multiplying the k and constant terms by L /R 2 and L^/R 4 , respectively. 

We can easily obtain an explicit expression for the wave functions themselves if we neglect 
the Pi, j s i and higher order terms, as the potential ( |5.30| ) is then that in the radial wave equation 
for a four-dimensional harmonic oscillator. To this order, the wave functions are given up to a 
normalization constant by 



^ = z 3 ' 2+l L^ +l) ( =-Qekz 



exp } — -Qekz 2 



(5.41) 



where, as before, v = (n — £)/2 is the order of the generalized Laguerre polynomial L^ +1 \ 

The dispersion relations (|5.38| ) are the central result of Section 5. We shall analyze (|5.38 ) and 
discuss its consequences at length in Sections 5.4 and 6. First, however, we close this more technical 
discussion with a few remarks related to the approximation that we have used to obtain the large-fc 
dispersion relations: 



1. The wave function is localized at the tip of the brane, near p = which is the fixed point of 
the 50(4) symmetry at which the S 3 shrinks to zero size and the fluctuations are fluctuations 
in R 4 . This is the reason why we find a four-dimensional harmonic oscillator. 

2. Our approximation is valid for wave functions that are tightly localized near z — 0. Evidently, 
this approximation must break down for mesons with high enough n, whose wave functions 
explore more of the potential. More precisely, if we increase n and £ while keeping v fixed and 

small, the wave functions are peaked at z ~ (^) 2 with a width 1 - T . Or, if we increase n 
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and v while keeping t fixed and small, the wave functions become wider, with v oscillations 

over a range of z from near zero to near z n ~ (ttt) 2 and hence a wavelength ~ — - — r . In 

either case, our approximation must break down for n ~ k, since for n this large z is no 
longer small and the wave function is no longer localized near z = 0. 

3. We must ask at what k (or, at what to) stringy effects that we have neglected throughout may 
become important in the dispersion relations for the mesons that we have analyzed. We can 
answer this question by comparing the length scale over which the meson wave functions that 
we have computed varies to the string length scale a' 2. Considering first the case where v is 
small, we see from (|4.18|) that the proper distance between the maximum of the wave function 



at z = zq and the tip of the brane at z = is 

'■~^~^f«(s)' (5 - 42) 



and the width of the wave function is 



61 ~ ) £ R (-^- ) . (5.43) 



Stringy effects can be neglected as long as 61 ^> a' 2, meaning 

k < O(X^M) , (5.44) 
where in the last expression we have restored the dimensions of k using ( p.20|) and (|3.34j ). 



(Since to = Vgk at at large /c, this parametric criterion is the same for u as for k.) If v is large, 
the wavelength of the wave function should be compared to ct 'a meaning that 61 is reduced 
by a factor ~ an d stringy effects can be neglected only as long as 

k < 0(\*M/v) . (5.45) 

We can conclude from either (|5.44|) or ( |5.45| ) that we are justified in using the dispersion 
relation that we have derived in the k — > 00 limit, as long as we take the A — > 00 limit first. 6 

4. Notice that as e — > 1 (i.e. approaching the critical embedding), both vq and Q vanish. Our 
approximation will therefore break down at the critical embedding. (One way to see this is to 
note that in the leading terms in ( |5.30| ) we will then have zero times infinity, meaning that it is 



no longer obvious that these are the leading terms.) However, the first order phase transition 
occurs at e = 0.756, long before this happens. 

5.3 Numerical results 

We can also obtain the meson wave functions and dispersion relations numerically, without making 
either a small e or a large-/c approximation. In this subsection we plot a few examples of such 

6 Recall that although the mesons that we have focussed on have masses of order M ~ m q j\f\ 1 there are also 
higher-lying stringy mesonic excitations with masses of order M\i ~ m q /\i. Requiring A 1 / 4 to be large is what 
justifies our neglect of these stringy mesonic excitations, just as it justifies our neglect of stringy corrections to the 
dispersion relations of the low-lying mesons. Note also that the latter becomes important at an u> of order the mass 
of the former. 
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Figure 7: Potential and ground state wave function for (left three panels) and ip2 (right three panels) 
for k given by 5, 20 and 100 (top to bottom). All plots have e = 0.756, corresponding to the Minkowski em- 
bedding at the dissociation transition. V{z) and the ground state (n = i = 0) solutions to the Schrodinger 
equation in the potentials V are both shown as solid lines, and the ground state energies are indicated by 
the horizontal (red) lines. The dashed lines show the approximation (|5.41| ) to the wave functions. 
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Figure 8: Dispersion relations for the ground state ipi meson with n 



at various values of e (i.e. 



771/ 



at various temperatures). The top (red) curve is the zero temperature dispersion relation uj 
with m given by ( [3.17 ) and with a group velocity that approaches 1 at large k, as required in vacuum 
by Lorentz invariance. The next three solid (black) curves are the dispersion relations for e = 0.25, 0.5 
and 0.756, top to bottom, the latter corresponding to the Minkowski embedding at the temperature T^ iss 
at which the first order phase transition occurs. The dashed (red) lines are the large-A: approximation 
discussed in Section |5.4| , given by uj(k) = v^k + Q,eLq/{vqR 2 ) with Q, specified by (|5.53| ). We see that the 
dispersion relations approach their large-A: linear behavior from below. The limiting velocity vq decreases 
with increasing temperature. Had we plotted dispersion relations for 0.756 < e < 1 corresponding to 
metastable Minkowski embeddings with T > Td; ss , we would have seen vq — > as e — > 1, approaching the 
critical embedding. 



results, and compare them to the analytic expressions that we have derived above upon making the 
large-fc approximation. 

In Fig. [7| we plot the potentials (|4.22|) and ground state wave functions for those potentials 
that we have obtained numerically for three values of k. Note the changing vertical scale in the 
plots of V; as k increases, V deepens. We see that as k increases and the potential deepens, the 
wave function gets more and more localized near z = and, correspondingly, the expression ( f5.41| ) 
for the wave function that we have derived in the large-/c limit using the fact that the wave function 
becomes localized becomes a better and better approximation to the exact wave function. 

In Fig. M we show dispersion relations obtained numerically for the ground state ip\ meson at 
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Figure 9: Group velocities v g = dto/dk for the dispersion relations from Fig. ^, with e = 0.25, 0.5 and 
0.756 (top to bottom). We see that the group velocity approaches its large-/c value vo from above. And, 
we see vq decreasing with increasing temperature. (Again, vq would approach zero if we included the 
metastable Minkowski embeddings with T > Tdi ss .) 

several values of the temperature. At each k, we solve the Schrodinger equation to find the ground 
state (using the shooting method) and from the eigenvalue we obtain u 2 and hence a point on the 
dispersion relation. By doing this at many fc's, we obtain the curves plotted. We also overlay the 
linear approximation to the large- k dispersion relations that we shall discuss in Sect. |5.4| . In Fig. H 
we plot the corresponding group velocities. 

5.4 Summary, limiting velocity and dissociation temperature 

In this Section we restate our central result for the dispersion relation and then discuss its implica- 
tions vis a vis a limiting velocity for mesons at a given temperature as well as a limiting temperature 
below which mesons with a given velocity are found, and above which they are not. 

In Section 5.2, we have derived the large-/c approximation to the meson dispersion relations at 
any temperature below the dissociation transition. We have checked this result against numerical 
solutions valid at any k in Section 5.3. We begin by restating the analytic result ( |5.38| ): 



co 2 = v 2 k 2 + ns(n + 2)^k + d snl % , +0(l/k) (5.46) 

K z it 4 



where 



i-s 2 Q2 32(l- g 2 ) 2 (l + 5 4 ) 

^ = T+F' n= (TT^ ' (5 ' 47) 
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The constant term d sn £ (which depends on whether we are discussing the ipi or ip 2 mesons — s — 1 
or s = 2 — and on the quantum numbers n and t) was given in ( |5.39| ) and ( |5.4(J| ). In writing 



the dispersion relation ( [5 . 46|) we have restored dimensions by undoing the rescaling (|3.34|) . The 
dimensionful quantity that we had scaled out and have now restored can be written as 



L ( 27im q \ [Too 



R 2 V V\ J V £ 




(5.4* 



where we have used ( |3.23| ), ( |3.32| ) and ( |3.33| ). The first factor in ( |5.48| ) is a (dimensionful) constant. 



The quantity too/e appearing in the second, dimensionless, factor is weakly temperature dependent: 
it can be read from Fig. ^ and is not constant to the degree that the curve in this plot is not a 
straight line (in the relevant regime < e < 0.756, as e = 0.756 corresponds to T = Tdi SS .) Although 
using dimensionless variables obtained via scaling by the temperature-dependent L /R 2 was very 
convenient in deriving all our results, in plotting the dispersion relation and group velocity in Figs. [8| 
and H we have instead plotted uj and k in units of T diss = 2.166 m q / y/X, which is a relevant, constant, 
physical, quantity comparable in magnitude to Lq/R 2 . In the remainder of this section, we shall 
analyze Q5.46j ). 

In the large-fc limit, the asymptotic value of the group velocity du/dk is given by vq. This 
velocity decreases with increasing temperature, and vanishes as e — > 1 on the critical embedding 
that separates Minkowski and black hole embeddings in Figs. [I] and [2| At the temperature at which 
the first order dissociation transition occurs, e = 0.756 and vq = 0.273. 

There is a natural explanation within the dual gravity theory for how the asymptotic velocity 
vq can arise. Using ( 3.37 ), it is easy to show that vq in (|5.47|) can also be written as 

which we see from ( |3.35| ) is precisely the local speed of light at the tip of the D7-brane. (The local 
speed of light is 1 at u = oo, and decreases with decreasing u, decreasing to vq at the tip of the 
D7-brane where p = and u — y = 1.) Since we have seen that in the large- A; limit the wave 
function of the meson fluctuations becomes more and more localized closer and closer to the tip 
of the D7-brane, this makes it natural that Vq emerges as the asymptotic velocity for mesons with 
large k. 

In the low temperature (equivalently, heavy quark) limit, we find (either directly from ( |5.47| ) 
or, initially, in (|5.27 ) in Section 5.1) that 



v 2 w 1 - As 2 . (5.50) 
Since ~ e at small e, using (|3.32j) we have 



2 A 2 T 4 , . 

v 2 ^l , 5.51 

16m 4 ' v ; 



which is precisely the critical velocity ( |2.7| ) obtained in [23| from the screening calculation as the 
velocity above which the potential between two moving quarks of mass m q cannot be defined. This 
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Figure 10: The asymptotic velocity vq from Q5.47 ) as a function of e. The low temperature approximation 
( 5,50| ) is plotted as a dashed line. Recall that the dissociation transition occurs at e = 0.756. 



is the first of two quantitative comparisons that we will be able to make between our present results 
for meson propagation and results obtained previously via the screening calculation. We see from 
Fig. that ( |5~50D works very well where T <C m q /^/\, which is where it was derived (both here 



and in @). 

In order to analyze ( |5.46|) beyond the k 2 term, it is instructive to rewrite it as a large-A; 
approximation to the dispersion relation u itself rather than to u 2 , yielding 

ne(n + 2)L Ad sn£ v 2 - n 2 e 2 (n + 2) 2 L 2 1 2 
= Vok + 2v R 2 + 8^ -W k + ° {1/k ] ' (5 ' 52) 



in the form we discussed in Section 1. We see that the term linear in k in (|5.46| ) yields a constant 



shift in the meson energies in ( |5.52| ). Whereas vq is independent of s, n and £, the constant term 
in ( |5.52|) results in evenly spaced dispersion relations for mesons with differing n quantum number, 



(leL /arg. ( ^ 3) 



separated by 

2r l) R- V vA 

which we plot in Fig. [11]. 

If we neglect the 0(1/ k) and higher order terms in (|5.52|) , the dispersion relations are the 




same for mesons ipi and ip2 and are independent of £. These degeneracies are broken at order 
0(l/k), where d sn £ first appears. We find that the coefficient of 1/k in uo(k) of (|5.52|) is typically 
negative: it is negative at all e < 1 if £ = for any n; it can become positive only if e, n and I 
are all large enough. When this coefficient is negative, it means that u(k) approaches its large-A; 
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Figure 11: The ^-independent spacing Q,eLq/2vqR 2 between the dispersion relations for any two mesons 
whose n quantum numbers differ by 1, in units of Tdi ss . See ( [5. 53 ). 



asymptotic behavior (which is a straight line with slope vq offset by the constant term in (|5.52|) ) 
from below. This means that d 2 uj/dk 2 < at large k and means that the group velocity v = duo/dk 
approaches vq from above at large k, as shown in Fig. [5[ However, at k = the group velocity 
vanishes and d 2 u/dk 2 > 0. (We have shown this analytically at small e in Section 5.1, see (|5.11| ), 
and our numerical results as in Section 5.3 indicate that this is so at all e.) So, as a function of 
increasing k, the group velocity begins at zero, increases to some maximum value that is greater 
than vq, and then decreases to v o as k — > oo as depicted in Fig. |9]. 7 Although vq is not the maximum 
possible group velocity, it appears that the maximal velocity exceeds vq only by a small margin. 
For example, for the ground state ipi meson whose dispersion relations are given in Figs. || and || 
we find that vq = 0.882, 0.6, and 0.273 for e = 0.25, 0.5, and 0.756 whereas the maximal velocities 
are 0.896, 0.634 and 0.342, respectively. We shall therefore simplify the following discussion by 
taking the maximal possible meson velocity at a given temperature to be the limiting velocity vq, 
neglecting the slight imprecision that this introduces. 

We now wish to compare our results for the limiting meson velocity vq at a given temperature 



to the result ( |2.6| ) inferred (qualitatively) from the analysis of screening in a hot wind in P3" |. We 
must first convert Vq(e) into Vq(T), meaning that we must convert from e to as discussed in and 
around Fig. |2|. The result is the solid curve in the left panel of Fig. [T2|, where we have plotted vq 
versus T/T^ ss . We have derived this curve as a limiting meson velocity at a given temperature. 



7 This behavior is not inconsistent with our identification of v with the local speed of light at the tip of the brane: 
it is only for k — > oo that the meson wave function is squeezed down to the tip of the brane; at finite k, the wave 
function is peaked where the local speed of light exceeds vq . 
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Figure 12: Left panel: The solid curve is the limiting velocity vq as a function of T/T^ ss , where Tdi ss is 
the temperature of the dissociation transition at zero velocity. The dissociation transition occurs at the 
dot, where vq ~ 0.273. The dashed curve is the approximation obtained by setting f(v) = 1 in ( |5.54| ). 
Right panel: f(v), the ratio of the solid and dashed curves in the left panel at a given v. We see that f(v) 
is within a few percent of 1 at all velocities. 



However, it can just as well be read (by asking where it cuts horizontal lines rather than vertical 
ones) as giving Tdi SS (v), the temperature below which mesons with a given velocity v are found and 
above which no mesons with that velocity exist. We see that Tdi SS (f) — > for v — > 1, the regime 
where Vq is given by (|5.51| ) and Td iss (v) is therefore given by ( |2.S| ). In order to compare our result 



for Td iss (f) at all velocities to (|2.6|) , we parametrize our result as 

T diss (v) = f(v)(l - *; 2 ) 1/4 T diss (0) . (5.54) 



In the left panel of Fig. |12] we compare our result (the solid curve) to (|5.54|) with f(v) set to 1, 



which is of course ( |2.(j| ). In the right panel, we plot f(v). We see that this function is close to 1 
at all velocities, varying between 1.021 at its maximum and 0.924 at v — 1. The weakness of the 
dependence of f(v) on v is a measure of the robustness with which the simple scaling ( |2.6| ) describes 
our result for the meson dissociation temperature at all velocities. 



6. Discussion and Open Questions 



We have used the AdS/CFT correspondence to compute the dispersion relation uj(k) for the heavy 
"quarkonium" mesons that exist as stable bound states in the strongly coupled plasma of M = 4 
SYM to which heavy fundamental quarks with mass m q have been added. In Section 4 we have 
introduced a new, and more geometrical, method of analyzing these mesons that has allowed us, in 
Section 5, to obtain the dispersion relations at large-/c analytically in the form ( 5.52 ), which we can 
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summarize as in Section 1 by writing 

u(k) = v k + a + j + ... . (6.1) 

We have computed a and b explicitly and analytically in Section 5, but at present we have no 
argument that the behavior of these coefficients, which depend on the meson quantum numbers, 
could teach us lessons that generalize beyond the particular theory in which we have computed 
them. On the other hand, the limiting large-/c meson velocity vq seems to encode much physics that 
may generalize to meson bound states in other strongly coupled gauge theory plasmas. 



Our explicit result is 



e 2 



v = -— -7T , (6.2) 



l + e 

, pir in Pi 



where e is related to = AT 2 / (8m 2 ,) as in Fig. ||. We see that v depends on the temperature 



(in the combination \/\T/m q ) but not on the meson quantum numbers. We see in Figs. |8] 
and H that v decreases with increasing temperature, becoming much less than 1 as the 
temperature approaches Tdi SS , the temperature at which mesons at rest dissociate. We see in 
these figures that the coefficient b in (|6.1|) can be negative, meaning that the group velocity 
approaches its large-A; value vq from above. Thus, vq is the limiting meson velocity at large 
k, but the maximal velocity occurs at finite k and is slightly larger than vq. We describe this 
quantitatively in Section 5, but it is a small effect and in this discussion we shall ignore the 
distinction between vq and the maximal velocity. 

We find that Vq, which in the gauge theory is the limiting velocity of the mesons that they 
attain at large k, also has a nice interpretation in the dual gravity theory. It is precisely 
the local velocity of light at the "tip" of the D7-brane, namely where the D7-brane reaches 
closest to the black hole. This is physically sensible because we have shown that the D7-brane 
fluctuations — i.e. the mesons in the dual gravity theory — are localized at the D7-brane tip 
in the large- k limit. 

At low temperatures or, equivalently, for heavy quarks we find 



v « 1 . (6.3) 

This is precisely, i.e. even including the numerical factor, the criterion for meson dissociation 
inferred from a completely different starting point in [p3fl . The logic there was that the 
screening length that characterizes the potential between a quark and antiquark moving with 
v > Vq is shorter than the quark Compton wavelength, meaning that if a quark and antiquark 
moving with v > Vq are separated by more than a Compton wavelength, to leading order in 
y/X they feel no attractive force. By inference, no mesons should exist with v > vq. We now 
see this result emerging by direct calculation of meson dispersion relations, rather than by 
inference. 
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We have a result for Vq(T), the limiting velocity beyond which there are no meson bound 
states, at all T < Tdi SS not just at low temperatures, see Fig. O. We can just as well read 



this as determining a temperature T^ ss (v) above which no meson bound states with velocity 
v exist. We find that up to few percent corrections, see Fig. 0, this is given by 

T diss (v) = (l-v 2 ) 1/4 T diss . (6.4) 

Once again, this is a result that was previously inferred from analysis of the velocity depen- 
dence of the screening length characterizing the potential between a quark and antiquark 
moving through the plasma [2IJ. We have now derived this result and the (few percent) cor- 
rections to it for the mesons whose dispersion relations we have explicitly constructed. We 
should also note that it is a slight abuse of terminology to call T diss (v) at v > a "disso- 
ciation" temperature: although it is a temperature above which no mesons with velocity v 
exist, if we imagine heating the plasma through Tdiss(^) we have not shown that any mesons 
present therein dissociate — they may simply slow down. The question of what happens in 
this hypothetical context is a dynamical one that cannot be answered just from the dispersion 
relations we have analyzed. 



As we discussed in Section 2, the result ( p.4| ) can be read as saying that no mesons with velocity 
v exist when the energy density of the strongly coupled plasma exceeds Pdiss(^) where, up to 
small corrections, 

Pdiss(t') = (1 - ^ 2 )Pdiss , (6.5) 

with pass the energy density at which mesons at rest dissociate. Correspondingly, the low 
temperature result (|6.3|) can be written as 

P 

± — Vq = constant , (6.6) 

Pdiss 



valid when p <C p^iss and Vq is close to 1. Thinking as in ||38|| , we can ask whether the 
same result holds in other theories. It will be interesting to address this question in (3 + 1)- 
dimensional gauge theories that are in various senses more QCD-like than M = 4 SYM. At 
present, however, we have only investigated the (p + l)-dimensional gauge theories described 
by N Dp-branes into which fundamental quarks, and hence mesons, have been introduced 



by embedding a Dg-brane [|56|, ^8], |29|. The Dp-branes fill coordinates 0, 1, . . . ,p. The Dg- 
brane fills the first d+ 1 of these coordinates 0,1, ... ,d, where d may be less than or equal to p, 
as well as q — d of the remaining 9 — p coordinates. In Appendix B, we sketch an investigation 
of those theories for which p — d + q — d = A. (The case that we have analyzed throughout 
the rest of this paper is p — d — 3, q — 7.) These theories are not conformal for p ^ 3, as 
their coupling constant A has dimension p — 3. It is convenient to introduce a dimensionless 
X eS = XT P ~ 3 . We have not repeated our entire construction for the Dp-Dg-brane theories. 
However, we expect that the wave functions for large-/c mesons will again be localized at the 
tip of the Dg-brane, and therefore expect that in these theories vq will again be given by the 
local velocity of light at this location. We compute this velocity in Appendix B. Assuming 
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that this is indeed the limiting meson velocity, we find 

l- £ (?-?)/ 2 

_ 1 + e (7-p)/2 

where e is given at small T/m q by 



1_ £ (7-p)/2^ 

t; o = ( ; I , 7 _ n s». ) > ( 6 - 7 ) 



cx ( JL ) = ± -£ . (6.8) 



m q j m q 

(Relating e to beyond the small T/m q limit requires solving the embedding equation given 
in Appendix B.) In these theories, the energy density of the plasma depends on parameters 
according to |55| 

p oc n 2 t p+1 \ (p ~ 3)/{5 ~ p) = jv 2 A (p-3)/(5 ~ p) T (14-2p)/(5 ~ p) , (6.9) 

and zero-velocity mesons dissociate at some energy density pdiss corresponding to £ — £diss 
where £diss = From these results we notice that at small e 

£ {7 ~ p)/2 oc — v oc , (6.10) 

m q F Pdiss 

meaning that the velocity vq of ( |6.7| ) can be written in the form ( |6.6|) for all values of p\ In 



Appendix B, we describe the verification that (|6.4j ) also holds, but only when phrased as in 



7SD in terms of energy density rather than temperature. 



Emboldened by these successes, we advocate investigating the consequences that follow from 
hypothesizing that T and J/^f mesons in the strongly coupled quark-gluon plasma of QCD propa- 
gate with a dispersion relation flB.ip with vq dropping dramatically as the temperature approaches 
Tdi S s from below, and with no bound states with velocity v possible if T > T diss (v) given by (KM ). In 
applying ( KM ) to QCD, it is important to scale T diss (t>) relative to the T d i SS for T and J/^f mesons 
in QCD itself. The result Tdi SS = 2.166m g /v^A for the mesons that we have analyzed is surely 
affected by the fact that they are deeply bound and so should not be used as a guide to quarkonia 
in QCD. For example, it seems to overestimate Tdi aa for J/ty mesons by a factor of 2 or 3. However, 
as argued in |22], |23| and as we have discussed above, the velocity scaling (|6.4|) may transcend the 
detailed meson physics in any one theory and apply to mesonic bound states in any strongly coupled 
plasma. The successful comparison of our detailed results to this simple scaling form supports this 
conjecture. 

As we have explained at length in Section 1, meson propagation is only one piece of the 
physics that must be treated in order to understand quarkonium suppression in heavy ion collisions. 
Introducing the dispersion relation and limiting velocity that we have found into such a treatment 
is something we leave to the future, instead making only a few qualitative remarks. 

First, from the dispersion relations alone we cannot conclude that if a quark-antiquark pair is 
produced (from an initial hard scattering) with a velocity v > Vq(T), with vq(T) the limiting meson 
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velocity in the plasma of temperature T in which the quark-antiquark pair finds itself, then the 
quark-antiquark pair do not bind into a meson. The reason that we cannot make this inference 
is that the dispersion relations describe stable mesons with arbitrarily large momentum k, making 
it a logical possibility that a high velocity quark-antiquark pair with arbitrarily high momentum 
interacts with the medium in some way such as to slow down and lose energy while conserving 
its momentum, and thus in some way dresses itself into a meson with arbitrarily high momentum 
k, and velocity v . That is, since the dispersion relations describe the propagation of mesons 
with arbitrarily large momentum, by themselves they do not require that quarkonium production 
is suppressed when the precursor quark-antiquark pair has velocity v > v (T). Excluding this 
possibility, allowed by the kinematics, requires some consideration of the dynamics. The heuristic 
argument of [^| provides guidance: the precursor quark-antiquark pair with v > Vq(T) do not 
attract each other and so even though it is kinematically allowed by the meson dispersion relations 
for them to slow down and form a meson, instead they will propagate independently through 
the medium. Thus, the pr-dependent quarkonium suppression pattern proposed in [p2| , with the 
production of quarkonium states with Tdi SS higher than the temperature reached in a given heavy 
ion collision experiment nevertheless becoming suppressed above a threshold transverse momentum 
at which a quark-antiquark pair with that transverse momentum has velocity Vq(T), rests upon 
the dynamical argument of ||23|| . It is natural that analyzing quarkonium suppression requires 
consideration of both the precursor quark-antiquark pair and the putative meson, and only the 
latter is described by the meson dispersion relation. It is then nice to discover that the limiting 
meson velocity vq(T) agrees precisely with the velocity at which quark-antiquark pairs can no longer 
feel a force at order \/X. 

We have just argued that the very large- /c region of the meson dispersion relation is unlikely 
to be populated in heavy ion collisions. But, whether or not such large-fc modes are excited, it is 
clear from Fig. [| that at temperatures near to Tdi ss mesons at any k move much more slowly than 
they would if they propagated with their vacuum dispersion relation. There are several in-principle- 
observable signatures of the slow velocity of quarkonium mesons. First, it increases the separation 
in space long after the collision between those mesons that are produced at the surface of the 
fireball moving outwards, and hence escape into vacuum promptly, and those which are produced 
in the center of the plasma and hence move more slowly than if they had their vacuum dispersion 
relation. An increase in the typical separation of identical mesons because of this slow velocity 
effect will shift the onset of Bose-Einstein enhancement in the two particle momentum correlation 
to a lower relative momentum. This simple idea underlies a technique widely used in heavy ion 
physics and often referred to as Hanbury-Brown Twiss (HBT) two-particle interferometry, in which 
identical two-particle momentum correlations are used to determine spatio-temporal characteristics 
of the collision region. For a review, see Ref. f57j. Quarkonium HBT interferometry would thus in 
principle be able to find signatures of a depressed meson velocity. Second, non-identical two-particle 
correlation functions are sensitive to whether one particle species A is emitted from the medium 
on average before or after another particle species B. Such a difference in average emission times 
could result, for instance, if the maximal velocities in the dispersion relations for A and B differ 
because of their different mass. The analysis of the effect of a difference in average emission times 



on non-identical two-particle correlation functions can be found in [58]. In principle, this provides 
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a second way of finding signatures of a depressed velocity for those mesons for which the plasma 
reaches temperatures close to their dissociation temperature. 

Quarkonium mesons in the quark-gluon plasma of QCD have nonzero width. In contrast, the 
mesons we have analyzed at T < Tdi SS are stable, with zero width. The dispersion relations that we 
have found have no imaginary part. This is certainly an artifact of the large number of colors N 
and large coupling A limits that we have taken throughout. Processes in which one meson decays 
into two mesons are suppressed by 1/N. And, thermal fluctuations which unbind a meson whose 
binding energy is 2m q are suppressed by the Boltzmann factor 

exp(-2m q /T) = exp(-0.92v / AT diss /T) , (6.11) 

which at some fixed T/T^ is nonperturbative in an expansion about infinite A. A calculation of the 
imaginary part of the meson dispersion relations at finite A remains for the future, but this simple 
consideration is enough to be sure that it is nonzero, as is the case in QCD at weak coupling [[II] . 



As soon as the mesons have nonzero width, their slow velocity has a further consequence in the 
context of heavy ion collisions: because they move more slowly, they spend a longer time in the 
medium giving the absorptive imaginary part more time to effect the dissociation of the meson than 
would otherwise be the case. 

Our discussion in this Section has highlighted three different avenues of further investigation 
opened up by our analysis of meson dispersion relations in a strongly coupled gauge theory plasma. 
The first is the investigation of the phenomenological consequences for J/^> and T suppression in 
heavy ion collisions of a dispersion relation of the form ( |6.1| ) with ( |6.4j ). Second, it appears to us that 
the most interesting open question about the mesons whose dispersion relations we have analyzed 
is extending the calculation to finite A and analyzing the width of the mesons. And, third, we could 
gain significant confidence in the application of the lessons we have learned to QCD by repeating 
our analysis for heavy quark mesons in the plasma of other strongly coupled gauge theories, in 
particular those with a controlled degree of nonconformality. 
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Appendices 

A. General discussion of brane embedding and fluctuations 

In this appendix we present a general discussion of brane embedding in a curved spacetime (in the 
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absence of fluxes) and its small fluctuations. We then specialize to the case of D7-branes embedded 
in the AdS$ x S5 black hole geometry. 



A.l General discussion 

Consider a p + 1-dimensional brane in a .D- dimensional target space whose action is 



Sd p = -Hp / <P +1 i J-dethij , (A.l) 



where £\i = 0,1, ... ,p denote the worldvolume coordinates and h^ is the induced metric in the 
worldvolume 

dX^ dX u 

kj = G^X) 9 d , n = 0, 1, . . . , D - 1 . (A.2) 



Suppose that Xq solves the equations of motion following from ( |A.1|) , thus describing an em- 



bedding of the brane in the target spacetime. We are interested in understanding the behavior of 
small fluctuations around Xq. For this purpose, let 

X»(t) = XS(?)+6X"(e). (A.3) 



The action for <5X M can then be obtained straightforwardly from ( A.l|) . The resulting action and 
equations of motion for 5X fl are, however, not geometrically transparent. This is due to the fact 
that SX^(^) is the difference between coordinates and thus does not have good properties under 
coordinate transformations. A more convenient way to parameterize 5X M (£) is to use the exponential 
map to express it in terms of a vector in the tangent space at Xq, as we now describe. (Such 
techniques have also been used in the calculation of string worldsheet beta functions Given a 



vector rf, we shoot out geodesies of unit affine parameter from Xq with tangent rf. The end point 
of such a geodesic is identified with Xq + 5X^. Such a map should be one-to-one within a small 
neighborhood of Xq. To second order in 77 one may solve the geodesic differential equation, finding 

6x» = rf- ^ryxo^V + • • ■ ■ (A.4) 

Note that the appearance of V is consistent with the coordinate dependence of SX; they can both 
be shown to have the same variation under a coordinate transformation. 



Using the parametrization (A.4), we find that 



hij = G^Xq + 8X)d i {X V Q + 5X v )dj(X* + 6X») (A.5) 

with 

h l3 = G^XoWiXgdjXZ = AfA J>; V, = \y„ Af = . (A.6) 
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The simplest way to find ( A. 5 ) is to use the Riemann Normal coordinates at X Q in which the 
Christoffel symbols vanish, hij is the induced metric in the worldvolume theory and below indices 
i,j will be raised and lowered by h. To quadratic order in rj we have 

-deth tJ = ^-deth tJ (l + XIV ^ + - (A^V^X^V'V) 

^W^ + ^^A^j • (A.7) 

We now take 77^ to be orthogonal to the brane worldvolume (which corresponds to choosing the 
static gauge), i.e. 

rf = XsnH{X ), s = l,...,D-p-l , (A.8) 

where n s (X ) are unit vectors orthogonal to the worldvolume direction. Note that Af and n% 
together span the full tangent space at X . i.e. 

Afn SM = 0, n s X = 5su KK = > (A.9) 

and 

^ = AfAt + <n s „. (A. 10) 

We now introduce 

K sij = Af \»V,n su , K s = K sij h ij , W st = <V^ = <AfV M r^ . (A.ll) 

K s ij is the extrinsic curvature of the brane in the s-direction, and is symmetric in (This follows 
from the fact that a surface orthogonal to satisfies V^n*, =^2 S v \u n l\ f° r some one-form u*. Note 
also that K si j can be written as K si j = ^L ns hij, where L n is the Lie derivative along n-direction.) 
U l st , which is antisymmetric in s,t, is an SO(D — 1 — p) connection for the transverse directions. 
Note that the choice of n£ (and thus Xs) is not unique. One can choose a different set of basis 
vectors by making an arbitrary local SO(D — 1—p) transformation. Thus Xs transforms as a vector 
under the SO(D —1—p) "gauge" symmetry and U l st transform as a connection. Note that this 
gauge symmetry is not dynamical. With these definitions we can now write 

Vify = (DiXs)n afl + KsijXstft , ( A - 12 ) 

where 

DiXs = diXs + UistXt (A. 13) 
is an SO(D — p — 1) covariant derivative. Using (|A.12|) in ( |A.7|) , we now find that 

Sd p = -V P J V-det/i^ (l + XsK s + ^D lXs D l X s + l - X sXt (-K sij K^ + R sijt h ij + K s K t )\ 

(A. 14) 
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with R S ijt = n^nfXfXjRa^p. For X to satisfy the equations of motion, the terms in (|A.14|) that 
are linear in the x's have to vanish. This implies that 

K s = K sij h ij = 0, s = 1, . . . , D - p - 1 . (A.15) 

These are the embedding equations for the background. Thus, the action ( |A.14j ) for the small 
fluctuations to quadratic order becomes 

Sop = fi P J (P +l i v^det^~ ^- l -D lX sD' l x s - l - X sXt {-K sij Kl j + R sijt 0)) . (A.16) 

We have used both the embedding equations ( |A.15| ) and the action for the small fluctuations ( |A.16| ) 
in Section 4. 

The action ([AlBD can be further simplified if satisfies additional constraints. For example, 
if is proportional to a Killing vector, then 

K sij = . (A. 17) 

This follows from the fact that satisfies V( M n^ = v^n^ for some v^. If in addition is a 
hypersurface orthogonal, i.e. if it satisfies V^n^i = w^n^ for some one form w^, then 

U ist = 0, for all t . (A. 18) 

We have used this simplification in Section 4. 

Finally, note that equation ( A.16Q was written using the coordinate split ( A.9 ). One can write 
it and other equations in a more covariant way by introducing 

V = h ij X^X ju , V = h ij X^X», hT = W^X) , (A.19) 

and using these objects in place of and Af in various places. h^ u = g^ v — n Sfl n su is the covariant 
induced metric on the brane and is the projector onto the worldvolume directions. 

A. 2 D7-branes in AdS 5 x S 5 black hole 

We now specialize to the case of D7-branes considered in the main text, where we have two transverse 
directions with 

where Ni t % are normalization factors. In this case U l st is proportional to the two-dimensional anti- 
symmetric tensor e st . It is easy to see that n v 2 is both hypersurface orthogonal and proportional to 
a Killing vector (since nothing depends on 0). We thus have K 2 ij = and U[ 2 = 0. The action 
( [A.16Q now reduces to the form we have used in Section 4, namely 

S D7 = fPr J d 8 ZV^h (l + ^{d^f + ^{dfoY + l -m\<j>\ + \m\<^ , (A.21) 
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where the "masses" are given by 



m-, 



m 



-Rn - R 2112 - KujKl 3 , 



(A.22) 
(A.23) 



— i?22 — -^2112 j 

with -R2112, and -R22 as defined in (|4.14j) . In writing ( |A.21| )- (|A.23| ) we have used the identities 

-Rst - Rent ~ Rs22U s, t = 1, 2 (A.24) 
and the fact that i?i2 = for the AdS*, x 5*5 black hole spacetime. We can also use the generalization 



R sij th ij = ny^Xt^R^p 



sec 



of the Gauss-Codazzi relation for a codimension two surface which we derive in Section [A. 3 
Eq. flOSD , to write 

K Uj Ki j = - &R + R- 2R u - 2i? 22 - 2R 2112 . (A.25) 



Therefore, ml in ( A.22|) can equivalently be written as 

m\ = R n + Rixv2 + 2R 22 + {S) R-R, 
which is the form that we used in Section 4. 

A. 3 Gauss-Codazzi relations for co-dimension 2 

Define the covariant derivative on the D7 brane as 



(A.26) 



(A.27) 



This is equivalent to the covariant derivative defined with respect to h^. We can now use D a to 
define the curvature of the D7-brane and then relate it to the curvature of the full space. Calculations 



similar to those in |60[ reveal that 

<% fc ' = P(R)J + (K s ) ik (K% - (K s ) jk (K s )\ 



(A.28) 



where s labels the two directions perpendicular to the brane and is summed over. P(R) is the 
projection of the full Riemann tensor onto the D7-brane, 



fj.ua 



(A.29) 



Taking further contractions of Eq. ( |A.28| ) with 5 1 a and h tk and using Eq.( |A.10| ) gives 

(8) i? = R- 2R SS - R tsst + K S K S - (K sij Ki j ) , (A.30) 

where s,t are both summed. In the case of interest, where K 2 ij = because n 2 is proportional to 
a Killing vector and where K s = is the embedding equation, we obtain ( A.25 ). 
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B. Dj9-Dg-Brane Theories 



It will be of interest in future to study the degree to which the meson dispersion relations that we 



have derived, together with their consequences like (j6\3| ) and (|6.4j) , change as one modifies the gauge 
theory to make it more QCD-like. In this appendix, we report on a check that we have mentioned in 
Section 6 in which the gauge theory is modified, albeit not in the direction of QCD. We consider the 
(p + l)-dimensional gauge theories described by N Dp-branes into which fundamental quarks, 
and hence mesons, have been introduced by embedding Nf Dg-branes 0, |48|, |49, 29| . The Dp- 



branes fill coordinates 0,1, ... ,p. The Dg-branes fill coordinates 0,1, ... ,d, where d < p, as well as 
q — d of the remaining 9 — p coordinates. In the large- N limit, the near horizon geometry of the 
Dp-branes is dual to a (p+ l)-dimensional supersymmetric Yang-Mills theory with 16 supercharges 
that is nonconformal for p ^ 3. We will restrict to p < 5. In the Nf/N — > approximation, the 
Dg-branes live in the background Dp-brane geometry, and their back-reaction on the geometry can 
be neglected. Strings which stretch between the Dg- and the Dp-branes are dual to Nf fundamental 
quarks in the gauge theory. We shall set Nf = 1. And, scalar mesons in the gauge theory are 
represented by fluctuations of the position of the Dq-brane. The specific case that we have analyzed 
throughout most of this paper is p = d = 3, q = 7. In this more general setting, as in the specific 
case, there is a dissociation transition at some T diss at which the spectrum of meson fluctuations 
changes from discrete to continuous. 

The background Dp-brane geometry is described by the metric |HH 



/ D \ (3-p)/2 / r (p-3)/2 \ 

ds 2 = R 2 (£) {-fde + r^/V p + r —^- (rfp 2 + dy- + P 2 dnl_ d _ x + y 2 dQl p _ q+d ) j 



(B.l) 



and the dilaton 



where 



L, 



(p-3)(7-p)/4 



^) g sr (P-W- P )/^ (B . 2) 



( V 7-P _ A7-p)/2) 2 

S = "- (7 -- )/2 ( " J , (B.3) 



u 7 ~p + e( 7 -p)/ 2 

(7-p)/2 _ ..-(7-p)/2 f-.r-p , _(7-p)/2 



r v t-n — u 



(J-p + £ V-p)/i) ; (B.4) 
u 2 = y 2 + p 2 , (B.5) 

and where we are using dimensionless coordinates as in ( |3.35| ). The black hole horizon is located at 
u = u = sfe. L specifies the position where the Dg-brane that we introduce will sit, as follows. 
We shall embed a Dg-brane described, in the absence of fluctuations, by a curve y(p) with the 
Dg-brane placed such that its tip is located at p = and y = Lq, and then use Lq to rescale metric 
coordinates such that the tip of the Dg-brane is at y(0) = 1. After this rescaling, the metric and 
dilaton are given by (|B.1|) and ( B.2 ). The holographic dictionary determines the coupling, number 
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of colors, and temperature in the gauge theory via 

A 
A 



rifiTT 3 ^- 3 )/ 2 



N 
T 



r(V) 

2 p - 1 n p - 2 g s a'^ 2 



(7_ p )2(5-p)/(7-p) (5 



47T 



(in 



p)/2fl-l / ^0 



(5-p)/2 



Note that A has dimension p — 3, making it useful to define the dimensionless coupling 

A fiff = XT p - 3 . 



(B.6) 
(B.7) 
(B.8) 

(B.9) 



The differential equation that specifies the shape of the embedding curve y(p) can be derived 
as we did in obtaining (|3.31| ). For the special case in which p — d + q — d = 4, the embedding 
equation simplifies, becoming 



y" 



+ 



(q-d- l)y' 2e^' 2 (y - y'p) ( (3 - d)u 7 -v + {q- d)e^ 7 



P 



-p)/2 



u 



2(7-p) _ e 7-p 



(B.10) 



We have scaled our variables so that the tip of the Dg-brane is at y(0) = 1; in order to have a 
smooth embedding we require y'(0) = 0; using these boundary conditions, we can then solve the 
embedding equation and obtain y(oo), which defines via y(oo) = y/e/e^,, ■ Finally, we can 
determine what the mass m q of the quarks that we are analyzing is via 



m„ 



eL 2 



47r 2 e 0O a' 2 



From (|B~6|) , @Bg ) and (|B~Tl| ) we find that 



-oo 



T 



A 



2/(5-p) 
off 



^2/(5-p) T 4/(5-p) 



(B.11) 



(B.12) 



where the constant a p is given by 

2 (10-2p)/(7-p) (3-p)/(5-p) 



(7 _ p)4/(5-p) 



r 



7 — p 



2/(5-p) 



We also note that the energy density of the plasma is given by |55] 



p _ & jV 2 T p+1 A^ -3)/(5-p) = 5 p A^ 2 A {p ~ 3)/(5 " p) T (14_2p)/(5 ~ p) 
where the constant b p is given by 



(9 - p) 2 6 -n^-^/^-p) 

(7_p)(19-3p)/(5-p) 



r 



7 — p 



2/(5-p) 



(B.13) 



(B.14) 



(B.15) 
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This means that 



-diss 



(7-p)/2 



P{T) 
Pdiss 



(B.16) 



where the zero-velocity mesons dissociate at a temperature Tdi SS corresponding to p = pdiss 

^.diss 

too 



and 



e£~, with 



,!ISS a constant of order unity. 

We shall not repeat our construction of the meson wave functions and dispersion relations for 
the Dp-Dg system here. Instead, we shall assume that in the large-fc limit the meson wave functions 
become localized at the tip of the Dg brane at p = and y = 1, as we found for the D3-D7 system. 
As a consequence, the limiting meson velocity will be given by the local speed of light at the tip of 
the Dg-brane. This velocity can be read from the metric ( |B.1|) , and is given by 



1 _ £ (7-p)/2 
1 + £ (7- P )/2 



(B.17) 



In Section 6 we have analyzed this result in the small e limit, showing that in this limit it takes 
on the form ( p.6| ) for any p. This illustrates the generality of the result ( |6.3| ) when it is phrased in 
terms of the energy density. Here, we shall analyze ([B.17| ) at arbitrary e < 1, seeking to compare 
it to (|6.4|). From ( |B.17| ) and ( |B.16| ) we see that the critical velocity satisfies 



1 - 1 -Vp 

1 + v ~ (1 + v ) 



e 



(7-J0/2 



P 



Pdis 



.diss 



(7-p)/2 



(B.18) 



Recall that e^, ss is a constant of order unity and that e/e^ is a, weak function of temperature and 
hence of p, obtained by solving the embedding equation and making a plot of 6^ vs. e as in Fig. ^ 
and reading off the ratio. 

Much as we did in Section 6, we can see (|B.18|) either as giving the limiting velocity Vq as a 
function of p, or as giving p<n S s{ v ), the energy density above which no mesons with velocity v exist, 
via 



Pdis 



[1-V )pdis 



'1 



.diss 



(P-7J/2 



(B.19) 



This is the generalization of ( p.4| ) to the Dp-Dg system. It is written somewhat implicitly, since e/e,^ 
which occurs within the square brackets is a weak function of Pdiss(f)- It is nevertheless manifest 
that the entire expression in the square brackets is a weak function of v, varying from one constant 
of order one at v = to some different constant of order one at v = 1. As in ( |5.54 ), we can then 
define a function f(v) by rewriting (|B.19|) as 



Pdi S 



r v \ = ry(^l(14-2p)/(5-p) Pdiss (0) 

7 2 



(B.20) 



where 7 = 1/ yl — v 2 is the Lorentz boost factor. Equivalently, using ( B.14| ) we can write 

Tdiss (0) 



T d i SS ( v ) = f( v ) 



7 (5-p)/(7-p) ' 



(B.21) 
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We have seen in Fig. 12 that for the D3-D7 brane system, f(v) is everywhere close to 1, with 
/(l) = 0.924 being the farthest it gets from 1. We have also done the exercise of solving the 
embedding equations for p = 4, the D4-D6 brane system with d = 3, and find in that case that the 
farthest that f(v) gets from f(v) = 1 is /(l) = 1.048. 

Given its derivation via QB.18 ), it would have been reasonable to try writing 



PdissW 



(14-2p)/(5-p) 1 - V 

Pdiss(0) TT^ 



(B.22) 



instead of ( |B.20| ). This does not work as well, yielding a f(y) that reaches 1.306 for the D3-D7 
system and 1.261 for the D4-D6 system. So although there is no important parametric difference 
between (|B.22j ) and ( |B.20| ), we have focussed on the form ( |B.2U| ), and hence ( |B.21| ), throughout this 
paper. 

The most important conclusion from our Dp-Dq investigation in this Appendix comes by com- 
paring ( p.20|) and ( |B.21|) . We see that in all the Dp-Dq systems we analyze, the leading velocity 
dependence of pdi SS (v) is that it is proportional to I/7 2 , as if the mesons see a boosted energy density 
as we discussed in Section 2. In contrast, T^ S8 (v) scales with a power of 7 that depends on p. 
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